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FOREWORD 


This  report  summarizes  the  comments  on  parachute  inflation  time  and  distance 
from  several  of  the  author's  previous  publications.  Additional  information  is  included 
to  explain  how  recovery  systems  opening  shock  forces  are  determined  not  only  by  the 
specified  operational  velocity,  altitude,  and  trajectory  angle,  but  are  also  determined 
by  the  initial  program  requirements  of  weight,  rate  of  descent,  equilibrium  altitude, 
choice  of  parachute  type,  and  canopy  airflow. 

Two  accepted  axioms  on  the  constancy  of  inflation  distance  and  the  effect  of 
inflation  time  on  the  magnitude  of  the  opening  shock  force  are  analyzed.  The  results 
verify  that  the  inflation  distance  axiom  is  correct,  and  that  the  inflation  time-opening 
shock  axiom  is  not  correct  although  it  appears  to  be. 


W.  W.  WASSMANN,  Head 
Underwater  Weapons  Division 
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INTRODUCTION 


One  important  aspect  of  parachute  technology  is  the  determination  of  the 
opening  shock  force  under  specified  operational  conditions.  Two  often  heard  common 
axioms,  concerning  the  parachute  deployment  process,  are: 

AXIOM  1:  The  product  of  the  system  velocity  at  suspension  line  stretch 

and  the  inflation  time  is  a  constant  inflation  distance. 

AXIOM  2:  Parachute  opening  shock  force  is  a  strong  function  of 

inflation  time. 

The  inflation  distance  of  axiom  one  is  presently  determined  experimentally 
from  field  test  data.  The  graphing  of  opening  shock  forces  versus  inflation  time 
certainly  presents  a  picture  that  strongly  supports  axiom  two.  What  is  needed  is  a 
theoretical  analysis  which  can  make  clear  the  relationship  of  the  two  axioms  to  the 
retardation  system  design  requirements,  operational  altitude  and  velocity,  trajectory 
deployment  angle,  type  of  parachute,  and  airflow  through  the  canopy. 

A  method  of  calculating  the  inflation  reference  time  and  distance  of  solid  cloth 
parachutes  was  developed  and  shows  that  a  parachute  system's  performance  is  not 
determined  solely  by  the  specified  operational  altitude  and  velocity.  It  is  more 
fundamental  than  that.  Rather  a  system's  performance  is  also  determined  when  the 
basic  parachute  system  requirements;  i.e.,  rate  of  descent  and  altitude,  system 
weight  to  be  recovered,  type  of  parachute,  and  canopy  cloth  are  selected. 

It  should  be  obvious  that  the  inflation  effects  depend  on  the  basic  system 
requirements  as  well  as  the  operational  requirements.  Sometimes  I  have  the  feeling 
that  designers  treat  the  basic  requirements  as  a  separate  and  independent  set  of 
calculations  from  the  deployment  effects  and  that  the  "system  performance"  is 
limited  to  application  of  the  operational  altitude,  velocity,  and  trajectory  angle. 

This  report  treats  the  development  and  examination  of  the  inflation  distance 
and  reference  time  of  parachutes.  It  is  a  composite  of  the  individual  observations 
from  the  several  references  with  additional  development  and  comments  to 
demonstrate  that  axiom  one  is  true  and  that  axiom  two  appears  to  be  true,  but  really 
is  not. 
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APPROACH 


The  inflation  of  any  parachute  is  a  very  complicated  process.  During  the 
inflation  time  period  the  mouth  area,  pressurized  canopy  surface  area,  and  parachute 
drag  area  are  varying.  The  instantaneous  rate  of  airflow  into  the  canopy  mouth  and 
out  through  the  transient  pressurized  canopy  area  are  changing.  Also,  the  velocity 
profile  during  inflation,  Reference  1, which  affects  the  canopy  rate  of  airflow  and  the 
inflation  reference  time,  is  dependent  upon  the  system  weight,  operational  altitude, 
parachute  size,  and  trajectory  angle,  Reference  2.  Under  certain  conditions  of  canopy 
airflow  and  altitude,  the  parachute  may  partially  inflate.  The  progression  of  the 
inflation  process  is  a  function  of: 

a.  Type  of  parachute 

b.  Weight  to  be  recovered,  including  the  weight  of  the  parachute 

c.  Parachute  drag  area 

d.  Canopy  rate  of  airflow 

e.  Operational  altitude 

f.  Instantaneous  trajectory  velocity 

The  calculation  of  the  inflation  reference  time  is  required  for  particular 
solutions  to  opening  shock  problems.  Methods  using  the  listed  criteria  in  a  basic  force 
analysis  are  developed  in  References  3, 4,  and  5.  Any  meaningful  analysis  of 
inflation  forces,  inflation  reference  time,  and  distance  must  include  the  stated  effects, 
as  well  as  produce  the  effects  that  have  been  observed  in  the  field. 


DEVELOPMENT  OF  THE  INFLATION  TIME  AND  DISTANCE  EQUATIONS 

The  nomenclature  of  the  various  parameters  to  be  discussed  is  presented  on 
Page  16  of  Appendix  A.  (See  Appendix  B  as  a  guide  for  the  use  of  Appendix  A.) 

The  referenced  opening  shock  force  studies  used  ratios  of  drag  area,  velocity  and 
time  to  determine  the  shock  factors,  and  the  time  of  occurrence  of  maximum  shock 
force  during  the  inflation  process.  The  ratio  concept  is  an  ideal  method  to  analyze  the 
effects  of  the  various  parameters  on  the  velocity  and  force  profiles  of  the  opening 
parachutes;  however,  a  means  of  calculating  the  inflation  reference  time  to  is 
required  before  specific  values  can  be  computed.  Methods  for  computing  the  varying 
mass  flow  into  the  inflating  canopy  mouth,  the  varying  mass  flow  out  through  the 
varying  inflated  canopy  surface  area,  and  the  volume  of  air,  Y0>  which  must  be 
collected  during  the  inflation  process  are  required. 
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The  mass  of  air  flowing  into  the  canopy  mouth,  m  inflow,  has  the  velocity 
reduced  which  raises  the  canopy  internal  pressure.  This  pressure  differential 
relative  to  the  canopy  cloth  causes  some  of  the  air  to  escape  through  the  canopy 
surface,  m  outflow,  from  the  confinement  of  the  canopy.  The  rate  of  outflow  is 
determined  by  the  weave  of  the  cloth  for  solid  cloth  canopies  or  the  ribbon  grid 
spacing  for  geometrically  porous  canopies  and  the  pressure  differential.  Excessive 
rates  of  airflow  result  in  partially  inflated  or  "stpiidded"  parachutes.  This  is  usually 
considered  as  a  situation  to  be  avoided.  A  meaningful  inflation  analysis  must 
account  for  the  possibility  of  squidding.  The  mass  of  air  retained  witnin  the  canopy, 
dm,  is  the  difference  between  tne  inflow  mass  and  the  outflow  mass. 

dm  =  m  inflow  -  m  outflow 

p-^=pVAM  -  pPAs  (1) 


The  solution  of  Equation  (1)  consists  of  three  steps.  The  first  step  is  to  find  a  method 
of  expressing  the  transient  mouth  area,  Am,  and  transient  pressurized  canopy  area, 
As,  as  a  function  of  the  deployment  time  ratio.  Sections  I  and  II  of  Appendix  A 
describe  a  method  where  infinite  mass  wind  tunnel  deployments  of  several  types  of 
parachutes  were  used  to  determine  a  dynamic  drag  area  ratio  for  use  in  opening 
shock  analysis.  A  conclusion  derived  from  the  dynamic  drag  area  ratio  was  that  the 
geometry  of  an  inflating  parachute  was  independent  of  altitude  and  velocity.  Berndt 
and  DeWeese,  Reference  6,  tested  finite  mass  man  carrying  parachutes  at  various 
altitudes  and  velocities  and  had  the  same  results.  Their  data,  Figure  1,  was  plotted 
as  the  ratio  of  canopy  projected  area,  to  surface  area  Sp/S0,  versus  t/tf.  Some  very 
important  conclusions  can  be  drawn: 

1.  The  geometry  of  an  inflating  parachute  is  independent  of  altitude,  velocity, 
and  system  mass. 

2.  Since  one  study  demonstrated  the  independence  of  drag  area  and  another 
demonstrated  the  independence  of  the  projected  area  to  surface  area  ratio,  ALL  of  the 
transient  geometric  shapes  of  inflating  parachutes  are  repeatable  and  are  also 
independent  of  altitude,  velocity,  and  system  mass  restraints. 

3.  Since  the  inflating  geometry  is  repeatable,  the  transient  volume  of  air  that 
is  collected  is  also  repeatable. 

Figure  2  represents  a  solid  cloth  type  parachute  canopy  at  some  instant  during 
inflation.  At  any  given  instant,  the  parachute  drag  area  is  proportional  to  the 
maximum  inflated  diameter.  Also,  tne  maximum  diameter  in  conjunction  with  the 
suspension  lines  determines  the  inflow  mouth  area  (A- A)  and  the  pressurized  canopy 
area  (B-B-B).  This  observation  provided  the  basis  for  the  following  assumptions.  The 
actual  canopy  shape  is  of  minor  importance. 

a.  The  ratio  of  the  instantaneous  mouth  inlet  area  to  the  steady-state  mouth 
area  is  in  the  same  ratio  as  the  instantaneous  drag  area.  The  mouth  area  is  not 
necessarily  circular,  but  is  usually  an  irregular  shape. 

_  CpS  (2) 

Amo  CqSq 
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Generalized 


FIGURE  1.  NORMALIZED  CANOPY  AREA  GROWTH  DURING  INFLATION  OF  28-FOOT  (D0) 

SOLID  FLAT  CIRCULAR  PARACHUTE 
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FIGURE  2.  PARTIALLY  INFLATED  PARACHUTE  CANOPY 


b.  The  ratio  of  the  instantaneous  pressurized  cloth  surface  area  to  the  canopy 
surface  area  is  in  the  same  ratio  as  the  instantaneous  drag  area. 


S  _  CpS 
Sq  C^Sq 


(3) 


c.  Since  the  suspension  lines  in  the  unpressurized  area  of  the  canopy  are 
straight,  a  pressure  differential  has  not  developed,  and  therefore,  the  net  airflow  in 
this  zone  is  zero. 


Figure  3  illustrates  a  typical  solid  cloth  parachute  wind  tunnel  infinite  mass  force¬ 
time  history  after  snatch.  In  infinite  mass  deployment,  the  maximum  size  and 
maximum  shock  force  occur  at  the  time  of  full  inflation,  tf.  However,  tf  is 
inappropriate  for  analysis  since  it  is  dependent  upon  the  applied  load,  structural 
strength,  and  materials  elasticity.  The  reference  time,  to,  where  the  parachute  has 
attained  its  steady-state  aerodynamic  size,  CdSo,  for  the  first  time,  is  used  as  the 
basis  for  performance  calculations. 


-  UNFOLDING  PHASE  - - ELASTIC  PHASE 


TIME  RATIO  t/t0 

FIGURE  3.  TYPICAL  INFINITE  MASS  FORCE-TIME  HISTORY  OF  A  SOLID  CLOTH 

PARACHUTE  IN  A  WIND  TUNNEL 
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Infinite  mass  opening-shock  signatures  of  several  types  of  parachutes  are 
presented  in  Figures  2  through  6  in  Appendix  A.  Analysis  of  these  signatures  using 
the  force  ratio,  F/Fs  versus  time  ratio,  Vto,  technique  indicated  a  similarity  in  the 
performance  of  the  various  solid  cloth  types  of  parachutes  which  were  examined.  The 
geometrically  porous  ring  slot  parachute  displayed  a  completely  different  signature, 
as  was  expected.  These  data  are  illustrated  in  Figure  4.  If  an  initial  boundary 
condition  of  CdSi/CdSo  =  0  at  time  t  =  0  is  assumed,  then,  the  data  can  be 
approximated  by  fitting  a  curve  of  the  form. 

CDS  l  t  \6 

C^=k)  ;T1=T  =  0  (4) 


FIGURE  4.  DYNAMIC  DRAG  AREA  RATIO  VERSUS  TIME  RATIO 


Reference  3  adapts  the  drag  area  ratio  to  a  wider  range  of  applications  by  not 
assigning  a  particular  value  to  the  exponent.  The  form  was  introduced  where 


CpS  _  /_t_\j 
CpS0  \ tQ/ 


;  r|  =  r  =  o 


(5) 
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j  =  1  is  indicative  of  geometrically  porous  parachutes  and  j  =  6  represents  the  solid 
cloth  family  of  parachutes. 


The  second  step  in  the  solution  of  Equation  (1)  is  to  develop  an  expression  for  the 
canopy  rate  of  airflow.  Parachute  canopy  airflow  is  discussed  in  another  section  of 
this  report.  The  rate  of  airflow,  P,  through  various  canopy  surface  configurations  can 
be  expressed  by  the  formula: 


P  =  k(AP)n 
AP  =  Cpq 


(6) 


then: 


aP  =  C 
P  =  k  | 


pT  PV 
CppV 
2 


2 

2 


(6a) 


A  basic  mass  llow  equation  for  the  solution  of  the  inflation  time  and  distance  is 
obtained  by  substituting  Equations  (2),  (3),  (5),  and  (6a)  into  Equation  (1). 


P 


dV 

dt 


= pVAmo 


(7) 


Measured  values  of  n  indicate  a  data  range  from  0.555  through  0.771.  A 
convenient  solution  to  the  reference  time  equation  for  solid  cloth  parachutes  evolves 
when  n  is  assigned  a  value  of  1/2,  and  j  =  6. 


Yo 

dY=  Amo 
o 


(8) 


The  third  step  in  the  solution  of  Equation  (1)  is  to  determine  the  volume  of  air  to 
be  collected  during  inflation.  The  volume  of  air.Yb,  associated  with  the  steady  state 
fully  inflated  parachute,  is  determined  using  the  methods  of  Section  VIII  of  Appendix 
A.  By  definition  to  is  the  particular  instant  during  the  inflation  process  when  tne 
transient  parachute  drag  area  is  equal  to  the  steady  state  drag  area  for  the  first  time. 
The  time  required  to  collect  the  Yo  volume  of  air  is  a  function  of  the  parachute  type 
and  geometry,  canopy  airflow  properties,  operational  altitude,  and  trajectory 
velocity.  The  trajectory  velocity  variation  during  inflation  is  influenced  by  the 
system  weight,  parachute  drag  area,  operational  altitude,  and  trajectory  angle  at  the 
line  stretch  velocity,  Vs.  Reference  4  applied  Equation  (8)  to  a  vertical,  toward  the 
earth  trajectory.  Inference  5  applied  Equation  (8)  to  a  general  trajectory.  Both  of 
these  references  supply  computer  programs  which  use  numerical  integration  to 
calculate  the  inflation  reference  time,  to,  and  the  opening  shock  force-deployment 
time  profile.  Reference  3  limited  the  trajectory  to  a  horizontal  direction.  This 
approach  develops  an  analysis  that  does  not  require  the  use  of  a  computer.  Solutions 
are  readily  obtained  with  a  desk  calculator  for  convenience.  The  horizontal  force 
analysis  derives  closed  form  equations  for  the  inflation  velocity  profile,  shock  factors 
and  time  of  occurrence  of  the  maximum  force.  The  force  analysis  of  Appendix  A 
yields  the  velocity  equation. 


V  = 


;ti  =  7  =  0 


(9) 
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With  the  insertion  of  Equation  (9)  into  Equation  (7)  for  j  =  6  yields. 


2n 

dt  (10) 


The  development  of  Equation  (10)  from  Equation  (1)  has  expressed  the 
transient  canopy  volumes  as  a  function  of  inflation  reference  time.  This  avoids  the 
need  to  model  or  calculate  specific  intermediate  geometric  shapes.  The  only 
geometric  shape  that  must  be  defined  is  the  steady-state  volume,  lo>  and  this  is 
available.  Equation  (10),  for  a  general  value  of  n,  requires  a  computer  solution.  Page 
13  of  Reference  4  defines  important  calculation  considerations  to  be  used  in  the 
solution. 


A  convenient  solution  to  the  reference  time  equation  evolves  when  n  is  assigned 
a  value  of  1/2.  Integrating  Equation  (10)  and  using : 


Vst0  M  = 


2W 


pgCDS0 


14W 


V sPQCpSo 


gpy? 

2W 


•i 


(ii) 


Multiplying  both  sides  of  Equation  (11)  by  Vs  demonstrates  that  Vsto  =  A  constant 
which  is  a  function  of  altitude 


v  t  14W 

gpio 

(12) 

2W 

^MO~^SO  k  | 

m 

i  J 

Vs'°-  pgCDs0 

-l 

Equation  (12)  verifies  axiom  1. 


As  a  demonstration  of  Equation  (11)  the  35-foot  D0  diameter,  30-gore  solid  cloth  cloth 
parachute  of  Example  1  of  Reference  4  is  to  be  applied.  Deployment  conditions  are: 


a.  System  weight  =  200  lb. 


b.  Density  altitude  =  3,000  ft. 


c.  Air  density  =  0.0021753  slugs/ft3. 

d.  Line  stretch  velocity  =  340  fps 


e.  Drag  coefficient  =  0.75 

f.  Initial  drag  area,  r  =  0. 


g.  Trajectory  angle,  y  =  0° 
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h.  Parachute  surface  area  S0  =  962.1  ft2 

i.  Parachute  drag  area  CdSo  =  721.6  ft2 

j.  From  Table  2,  Page  A-15  of  Appendix  A,  the  following  inflated  shape  data 
were  obtained  for  a  30-gore  solid  cloth  flat  circular  parachute.  See  Figure  24  of 
Appendix  A  for  shape  nomenclature. 

,o  _  _ 14x200 _ |eK1_il 

°  0.0021753x32.2x340x721.58  L  J 

k.  Steady-state  inflated  canopy  radius,  a  =  11.69  ft 

l.  Steady-state  mouth  area,  Amo  =  399.53  ft2 

m.  Steady-state  canopy  volume  of  air  to  be  collected,  Yo  geometric  = 
4690.83  ft3 

The  MIL-C-7020,  type  m  canopy  cloth  airflow  constant  is  k  =  1.46  and  the 
average  pressure  coefficient  is  taken  to  be  Cp  =  1.7  then: 


K1  = 


32.2  x  0.0021753  x  4690.83 


2  x  200 


tQ  =  0.773  sec. 


721.58 


399.53-962.1  (1.46) 


1 0.0021 753  x  1.7  j1/2 


(ID 


Knacke,  Reference  7,  presents  the  inflation  time  as 
nDo 

'"vr 


8x35 
"  340 


(13) 

where  n  =  8  for  solid 
cloth  parachutes 


ti  =  0.824  sec. 

While  Equation  (13)  is  more  convenient  than  (11)  to  use,  Equation  (11)  gives  the 
effects  of  altitude,  mass,  and  cloth  rate  of  airflow;  both  values  are  estimates  of  the 
inflation  time.  Equation  (11)  was  theoretically  derived.  Equation  (13)  was 
empirically  developed  by  other  experimenters  from  numerous  field  test  data.  The 
parachute  system  of  Example  1  deployed  at  an  altitude  of  1000  feet  yields  values  of 
to  =  0.839  second  and  tf =  0.824  second.  The  closeness  of  the  theoretical  to  and  the 
empirical  tf  values  indicates  that  the  application  of  the  criterion  listed  below  is  a 
reasonable  approach  to  the  inflation  time  calculation. 


a.  Average  pressure  coefficient 

b.  Cloth  permeability 


c.  The  trajectory  velocity  is  taken  as  the  canopy  inflow  velocity. 
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d.  The  ratio  of  the  instantaneous  mouth  area,  Am,  to  the  steady-state  mouth 
area,  Amo.  is  equal  to  the  instantaneous  drag-area  ratio. 

e.  The  ratio  of  the  instantaneous  pressurized  canopy  area,  S,  to  the  canopy 
surface  area,  S0,  is  equal  to  the  instantaneous  drag-area  ratio. 

f.  The  volume  of  air  to  be  collected  during  the  inflation  process  is  represented  by 

lb. 


EFFECT  OF  DESIGN  REQUIREMENTS  ON  PERFORMANCE 


At  the  beginning  of  a  parachute  recovery  system  design,  basic  design 
requirements  are  laid  down  that  determine  the  aerodynamic  and  geometric  size  of 
the  decelerator.  These  requirements,  together  with  decisions  on  the  type  of 
parachute  to  be  applied,  and  canopy  rate  of  airflow  influence  the  operational 
performance  of  the  design.  Usually  the  basic  requirements  include  the  desired  rate  of 
descent,  Ve,  the  altitude  at  which  it  is  to  be  achieved,  p,  and  the  payload  weight,  W. 


Once  these  decisions  are  firm  the  system  W/CdSo  has  been  established.  At  the 
equilibrium  velocity  the  parachute  drag  force  is  equal  to  the  system  weight. 


W  =  1/2pVe2  Cd  S 
W  2 

5^7  =  1/2pv* 


o 


(14) 


The  particular  system  weight  to  be  recovered  defines  the  parachute 
aerodynamic  size  (CdSo)-  The  geometric  size  (Do)  and  inflated  shape  and  volume  are 
a  result  of  the  type  of  parachute  and  canopy  rate  of  airflow  selected. 

„  0  2W 

CDSo  =  — -  (15) 

pVe 


The  inflation  distance  is  dependent  on  the  parachute  geometry,  cloth  airflow 
properties,  W/CdSq.  and  altitude. 


Vsto  = 


14W 

pgcDs0 


(12) 


When  the  parachute  type,  canopy  cloth,  and  suspension  line  length  have  been 
selected,  the  only  variable  undefined  is  the  operational  altitude.  The  inflation 
distance  does  not  depend  on  the  velocity  at  suspension  line  stretch,  but  the  inflation 
reference  time  does.  Rather,  the  inflation  distance  has  been  preordained  by  the 
system  requirements,  and  initial  design  decisions.  For  a  given  recovery  system  and 
operational  deployment  altitude  the  inflation  distance  is  a  constant. 


Now  let's  examine  axiom  2.  Table  1  of  Reference  3  summarizes  the  parachute 

Serformance  formulas  for  several  values  of  j.  Parachute  performance  for  horizontal 
eployment  depends  on  the  type  of  parachute,  j,  initial  drag  area,  CdSi,  and  the 
Ballistic  Mass  Ratio  (BMR),  M,  which  is  not  a  function  of  the  parachute  velocity  at 
suspension  line  stretch,  but  does  depend  on  the  inflation  distance  and  operational 
altitude. 
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M  = 


2W 

pgvs  to  Qj  So 


(16) 


The  time  of  the  occurrence  of  the  maximum  shock  factor  during  the  finite  mass 
process,  and  the  maximum  shock  factor  are  defined  by  the  BMR.  For  solid  cloth 
parachutes;  j  =  6, let  i=  0. 


(17) 


xi  max 


(18) 


The  BMR  is  the  scale  parameter  which  determines  what  percentage  of  the 
steady  state  drag  force,  Fs,  is  to  be  realized  as  opening  shock  force. 


F  s  =1/2p\|cbS0 


(19) 


Example  1:  A  Test  of  axiom  2 

A  given  finite  mass  parachute  system  is  to  be  tested  at  several  operational 
velocities  at  a  given  altitude.  Examine  the  opening  shock  performance  and  graph  the 
variation  of  the  maximum  shock  force,  Fmax>  versus  inflation  reference  time,  to. 

What  conclusions  are  evident? 

Solution:  The  system  is  tested  at  an  initial  line  stretch  velocity,  Vs.  The 
inflation  distance  Vsto  from  Equation  (12),  is  independent  of  velocity  and  is 
determined  by  the  parachute  geometry,  altitude,  and  cloth  airflow.  The  inflation 
distance  is  used  to  determine  tne  BMR  of  Equation  (16).  The  time  of  occurrence 
during  the  inflation  process,  Equation  (17),  and  the  maximum  shock  factor, 

(Equation  (18)),  are  determined  by  the  BMR.  The  line  stretch  velocity  only  affects 
operational  performance  in  the  steady  state  drag  calculation,  (Equation  (19)).  The 
maximum  shock  force,  Fmax  is  given  by: 

^max  ~  xj  Fg  (20) 

max 

For  a  given  altitude  and  line  stretch  velocity  Vs  the  inflation  distance  is 
determined  by  Equation  (12)  and  Fmax  by  Equation  (20).  If  the  test  velocity  is 
doubled  the  inflation  distance  is  unchanged.  This  means  that  the  BMR,  time  of 
occurrence,  and  maximum  shock  factor  remain  unchanged.  However,  the  steady- 
state  drag  has  increased  by  a  factor  of  four.  The  maximum  shock  force  increases 
because  the  constant  shock:  factor  of  the  system  is  claiming  the  same  percentage  of  a 
much  higher  force.  From  Equation  (12)  tne  doubling  of  the  velocity  reduces  the 
inflation  reference  time  by  1/2.  If  the  test  velocity  is  reduced  by  1/2  the  inflation 
reference  time  increases  by  a  factor  of  two  and  the  maximum  shock  force  is  reduced 
by  a  factor  of  1/4  due  to  the  reduction  of  the  steady-state  drag  force. 

These  data  are  plotted  in  Figure  5.  First  examination  of  Figure  5  would  seem  to 
demonstrate  that  the  maximum  force  is  a  strong  function  of  inflation  time.  However, 
the  maximum  shock  force  varied  because  the  steady  state  drag  force  varied  due  to 
changing  Vs,  not  because  the  inflation  reference  time  was  altered. 
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FIGURE  5.  APPARENT  VARIATION  OF  MAXIMUM  PARACHUTE  OPENING  SHOCK  FORCE 

WITH  INFLATION  REFERENCE  TIME 


The  conclusion  to  be  drawn  from  Figure  5  is  that  axiom  2  appears  to  be  true,  but 
is  not.  The  inflation  distance  is  more  important  to  the  inflation  process  than  the 
inflation  reference  time. 

Since  the  inflation  distance  is  constant  V8  can  be  expressed  as: 

Constant 


in  which  case  Equation  (20)  can  be  expressed  as 

Constant  \2  (20a) 

/ 

This  approach  does  show  the  inverse  parabolic  variation  of  opening  shock  force 
with  inflation  reference  time,  but  to  is  a  dependent  variable  of  V8  which  is  not  a 
reversible  process.  So,  even  though  Equation  (20a)  is  mathematically  valid  it  does 
not  express  the  problem  in  a  usable  form.  It  is  physically  possible  to  specify  a  line 


F  —  x 

'max  Ai 


pCp  s0 


max 
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stretch  velocity.  The  specified  line  stretch  velocity  results  in  a  specific  inflation 
reference  time.  It  is  NOT  physically  possible  to  specify  only  an  inflation  reference 
time  and  have  this  result  in  a  specific  line  stretch  velocity. 


DISCUSSION  OF  THE  VARIABLES  OF  THE  INFLATION  REFERENCE  TIME 
AND  DISTANCE  EQUATIONS 


gp^o 

o 

(/) 

Q 

o 

i _ 

V  t  14W 

2W 

A 

A  -A  ,M12 
mmo  mso  K 1  2 

pgcDsQ 

—1 

The  inflation  distance  (Vsto)  appears  to  be  a  velocity-time  effect.  Equation  (12) 
indicates  that  the  inflation  distance  actually  depends  on  system  mass  ratios,  area 
ratios  and  canopy  airflow  properties.  The  several  terms  of  Equation  (12)  each  have  a 
physical  significance  on  parachute  performance.  The  following  paragraphs  identify 
and  explain  the  meaning  of  each  term. 


Weight-to-drag  area  ratio.  W/CnSo 

The  identical  parachute  used  in  finite  mass,  intermediate  mass,  or  infinite  mass 
systems  experiences  a  variation  in  inflation  time  due  to  its  mode  of  use.  Low  weight- 
to-drag  area  ratios  in  the  range  of  0.3  suggest  low  rate  of  descent  finite  mass 
deployments  as  in  the  recovery  of  personnel.  High  rate  of  descent  infinite  mass 
applications,  as  in  ordnance  usage,  may  have  weight-to-drag  area  ratios  in  excess  of 
thirty  four.  Equations  (11)  and  (12)  are  applicable  to  all  modes  of  operation.  At  the 
time  to  all  inflations  have  collected  the  same  volume  of  air.  Finite  mass  inflations 
experience  a  significant  velocity  reduction  during  the  inflation  process  which  limits 
the  rate  of  airflow  into  the  canopy  and  extends  the  inflation  time.  Infinite  mass 
deployments,  on  the  other  hand,  can  be  nearly  constant  velocity  which  maintains  a 
higher  rate  of  mass  flow  into  the  canopy  and  decreases  the  inflation  time.  Actual 
constant  velocity  infinite  mass  operation  produces  the  minimum  inflation  time.  The 
inflation  reference  times,  for  the  parachute  of  Example  1,  are  plotted  in  Figure  6  as  a 
function  of  weight-to-drag-area  ratio  for  Equation  (11).  The  finite  mass  and  infinite 
mass  limiting  Ballistic  Mass  Ratios  for  solid  cloth  parachutes  are  referenced  on  the 
graph.  An  arbitrary  Ballistic  Mass  Ratio  value  of  10  was  selected  as  the  boundary 
between  intermediate  mass  and  infinite  mass  operation  for  solid  cloth  parachutes. 
Figure  9  of  Appendix  A  shows  a  minimal  velocity  reduction  at  t/to  =  1  for  a  BMR  of  10 
when  compared  to  true  infinite  mass. 

The  trends  of  Figure  6  show  that  the  principal  variation  of  the  inflation 
reference  time  occurs  in  the  finite  mass  mode  of  operation.  By  the  transition  from 
finite  mass  to  low  intermediate  mass,  a  near  minimum  inflation  reference  time  is 
reached.  This  quickly  converges  to  an  essentially  constant  value  for  the  high  end  of 
the  intermediate  mass  and  infinite  mass  regimes.  Similar  trends  for  other  types  of 
parachutes  are  a  reasonable  extension  of  the  premise  since  it  is  based  on  variable 
rates  of  flow  filling  a  constant  volume. 
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FIGURE  6.  VARIATION  OF  THE  INFLATION  REFERENCE  TIME  AS  A  FUNCTION  OF  SYSTEM  WEIGHT 
FOR  THE  FLAT  CIRCULAR  SOLID  CLOTH  PARACHUTE  OF  EXAMPLE  ONE 
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Coefficient  14W/(  Pg  CnSJ 

This  coefficient  is  the  ratio  of  the  system  mass  (W/g)  to  an  atmospheric  mass 
(pCdSo)  per  foot  of  length  of  the  Ballistic  Mass  Ratio  atmosphere  cylinder  of  Figure  8 
of  Appendix  A. 


Exponent  coefficient  gPYo/2W 

This  dimensionless  coefficient  in  the  exponent  is  the  ratio  of  the  mass  of 
atmosphere  associated  with  the  steady-state  inflated  canopy  (p  \o)  to  the  system  mass 
(W/g). 

Inflation  Time  and  Distance  Exponent  Dimensionless  Area  Ratio  Term 


The  dimensionless  area  ratio  expression  is  the  ratio  of  the  parachute 
aerodynamic  size  (CdS0)  to  the  effective  mouth  area  expressed  in  the  denominator. 
As  the  rate  of  flow  through  the  canopy  varies  with  k,  p,n,  or  pressure  coefficient,  the 
inflation  characteristics  are  modified. 


when: 


0  = 


AMo 


(21) 


A  critical  condition  exists  and  the  parachute  fails  to  fully  inflate.  As  the  effective 
mouth  area  approaches  zero,  to  increases  exponentially.  It  is  possible  for  the  inflation 
time  to  be  extended  to  the  point  that  the  parachute  system  is  not  practical  even 
though  the  effective  mouth  area  is  not  zero.  If  the  critical  parachute  system  is 
deployed  at  a  higher  altitude  the  air  density  reduction  causes  the  effective  mouth 
area  to  become  positive  and  the  canopy  inflates.  Imporous  canopies  (k  =  0)  always 
inflate.  These  phenomena  agree  with  field  test  observations. 


As  a  demonstration  of  the  cloth  rate  of  airflow  effects  in  Equation  (11), 
replacement  of  the  canopy  cloth  of  Example  1  with  a  high  rate  of  airflow,  three 
momme  silk  cloth  (k  =  7.43),  yields  an  in 
the  Example  1  operational  conditions.  T 
feet  and  the  system  would  appear  to  an  o 


- q -  j —  — 

lation  reference  time  of  to  =  99.76  sec.  under 
lis  time  exceeds  the  time  of  fall  from  3,000 
jserver  as  a  partial  inflation.  The  same 


system  deployed  horizontally  at  20,000  feet,  at  constant  dynamic  pressure,  has  an 
inflation  reference  time  of  to  =  1.5  sec. 


The  effective  mouth  area  is  a  clue  in  the  solution  of  the  parachute  critical 
velocity  phenomenon.  Another  part  of  the  problem  is  the  effect  of  suspension  line 
length  and  number  of  gores  in  the  canopy.  Consider  the  parachute  as  consisting  of 
two  main  elements.  A  canopy  which  generates  the  aerodynamic  force,  and  a  number 
of  suspension  lines  that  transmit  the  force  to  the  payload.  The  two  elements  are 
joined  at  the  canopy  skirt  hem  to  complete  the  parachute  assembly.  The  parachute 
inflation  characteristics  depend  on  the  force  distribution  at  the  canopy  skirt  hem,  see 
Figure  7.  Due  to  the  cone  angle  of  the  suspension  lines  there  is  a  radial  component  of 
the  suspension  line  force,  Fsl,  normal  to  the  canopy  centerline  that  tends  to  collapse 
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FIGURE  7.  CANOPY  SKIRT  HEM  GEOMETRY  ANI)  FORCE  DISTRIBUTION  OF  FLAT 
PARACHUTES  IN  FULL  STEADY  STATE  INFLATION 


the  canopy.  The  radial  component  of  the  canopy  aerodynamic  force,  Fca>  tends  to 
inflate  the  canopy.  At  full  or  partial  steady-state  inflation,  the  outward  radial 
component  of  the  canopy  inflating  aerodynamic  force  is  in  balance  with  the 
suspension  line  collapsing  normal  force  component.  Partial  inflations  produce  a 
pseudo  hem  where  the  inflated  and  uninflated  canopy  sections  meet.  During  the 
inflation  process  the  velocity  head  of  the  fluid  entering  the  canopy  is  reduced, 
resulting  in  an  increase  in  canopy  internal  pressure.  As  long  as  the  fluid  inflow 
exceeds  the  outflow,  additional  fluid  volume  is  added  to  the  canopy  and  the  internal 

Pressure  wave  advances.  When  the  outflow  is  equal  to  the  inflow,  the  velocity  of  the 
owing  fluid  is  not  reduced  and  additional  higher  pressure  fluid  is  not  added  to  the 
volume.  The  internal  pressure  is  sustained  at  the  partially  inflated  level. 

The  suspension  line  length  and  the  number  of  canopy  gores  contribute  to  the 
equilibrium  diameter  size  of  the  partially  inflated  parachute.  Lengthening  the 
suspension  lines  improves  inflation  characteristics  by  reducing  the  radial  force 
component  tending  to  collapse  the  canopv.  Shortening  the  suspension  lines  has  the 
opposite  effect.  Increasing  the  number  oi  gores  in  the  canopy  reduces  the  angle  <p, 
Figure  7,  between  the  tangent  force  in  the  cloth,  Ft,  in  the  main  seam  and  the  canopy 
mouth  diameter.  This  angular  reduction  uses  the  canopy  aerodynamic  force  more 
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efficiently  and  contributes  to  inflation.  In  the  limit,  angle  </>  approaches  zero  decree 
as  the  number  of  gores  approaches  infinity.  There  are  three  combined  contributing 
factors  to  parachute  critical  velocity. 

a.  Canopy  rate  of  airflow  -  Determines  the  ability  of  the  canopy  to  entrap  the 
volume  of  fluid,  Y0,  required  for  full  inflation.  See  the  cloth  rate  of  airflow  section  for 
the  effects  of  altitude,  k,  and  n. 

b.  Length  of  the  suspension  lines. 

c.  Number  of  gores  in  the  canopy. 

The  interrelationship  of  these  three  factors  control  the  inflation  characteristics  of 
each  design. 


EFFECTS  OF  ALTITUDE  ON  THE  INFLATION  REFERENCE  TIME  AND 
DISTANCE 

Numerous  simplified  filling  time  equations  have  been  published  in  many 
reports  and  will  not  be  repeated  here.  Generally  there  are  two  forms: 

*  nDo 

tf  =  (23) 

t,  -  (24) 

Figure  8  shows  that  when  Vs  is  raised  to  the  0.9  power,  the  result  is  nearly  the 
same  as  dividing  Vs  by  2.  So  Equation  (23)  can  be  simplified  by  multiplying  n  by  2 
and  setting  the  exponent  of  Vs  to  unity. 


SUSPENSION  LINE  STRETCH  VELOCITY.  EPS 

FIGURE  8.  EFFECT  OF  THE  FILLING  TIME  EQUATION  EXPONENT  0.9  ON  THE 
TRAJECTORY  VELOCITY  AT  PARACHUTE  SUSPENSION  LINE  STRETCH 
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The  accepted  values  of  n  for  Equation  (24)  are  8  for  solid  cloth  parachutes  and 
14  for  geometrically  porous.  Table  1  lists  canopy  fill  factors  for  various  types  of 
parachutes. 


TABLE  1.  CANOPY  FACTOR,  n,  FOR  VARIOUS  PARACHUTE  TYPES 


PARACHUTE  TYPE 

CANOPY  FILL  FACTOR,  n 

REEFED 

OPENING 

DISREEF 

OPENING 

UNREEFED 

OPENING 

SOLID  FLAT  CIRCULAR 

ID^ 

ID 

8 

EXTENDED  SKIRT,  10% 

16*18 

4-5 

10 

EXTENDED  SKIRT,  FULL 

16-18 

7 

12 

CROSS 

ID 

ID 

11.7 

RIBBON 

10 

6 

14 

RINGSLOT 

ID 

ID 

14 

RINGSAIL 

7-8 

2 

7 

RIBLESS  GUIDE  SURFACE 

4-6 

J/|D  =  INSUFFICIENT  DATA  AVAILABLE  FOR  MEANINGFUL  EVALUATION 
Reproduced  from  Reference  7 


Figure  9  compares  the  inflation  times  calculated  by  Equations  (11)  and  (24)  for  a  35- 
foot  Do,  T-10  type  canopy  deployed  at  a  sea  level  line  stretch  velocity  of  Vs  =  211  fps 
with  a  system  weight  of  200  lbs.  Figure  9  shows  the  time  variation  for  altitudes  up  to 
50,000  met  when  testing  at  constant  dynamic  pressure. 
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W  =  200  LB,  D0  =  35  FT;  C0S0  =  721.58  FT2; 
VSL  =  211  FPS;  V0  =  4690.8  FT3;  n  »  10 
Vs  =  vsl  VPsl/P 


'■o 

SEC 


1.42 
1  08 
0.84 
0.55 
0  38 
0  26 


FIGURE  9.  VARIATION  OF  INFLATION  TIME  FOR  A  T-10  TYPE  PARACHUTE  SYSTEM  AS  A 
FUNCTION  OF  EMPERICAL  AND  DERIVED  METHODS  OF  CALCULATION  FOR  VARIOUS 

ALTITUDES 

The  inflation  time  of  Equation  (24)  produces  an  apparent  reduction  in  inflation 
time  due  to  the  increase  in  Vs  to  maintain  a  constant  dynamic  pressure. 
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The  decrease  in  inflation  reference  time  in  Equation  (11)  is  due  to  the  decrease 
in  air  density.  For  constant  dynamic  pressure  as  altitude  increases  the  test  velocity 
must  be  varied. 


Substituting  for  Vs  into  Equation  (11)  and  simplifying 


At  all  altitudes  the  Vsl  remains  constant,  and  the  primary  variable  is  density. 
Other  possible  variations  may  be  due  to  drag  coefficient  rise  caused  by  a  decrease  in 
canopy  cloth  airflow  properties  and  possible  changes  in  geometry  caused  by  the  rise 
in  drag  coefficient.  Otherwise,  the  reference  time  is  only  a  function  of  system 
constants  and  altitude. 

As  actual  values  of  the  rate  of  fluid  flow  exponent  n  approach  1.0  the  inflation 
reference  time  and  distance  are  extended.  Equation  (7)  was  programmed  in 
Reference  5  for  values  of  n  other  than  0.5  and  the  effects  of  an  arbitrary  trajectory 
angle  at  deployment  for  j  =  6.  Figure  10  illustrates  the  effects  of  altitude,  n,  and 
trajectory  deployment  angle  for  a  constant  velocity  at  all  altitudes,  and  Figure  1 1 
illustrates  the  same  variables  for  a  constant  dynamic  pressure  at  all  altitudes. 
Figure  12  presents  the  effects  of  n  and  altitude  on  the  inflation  reference  time  for 
constant  velocity  and  constant  dynamic  pressure  for  horizontal  deployment.  At  sea 
level  the  extension  of  the  inflation  reference  time  due  to  a  realistic  n  =  0.632  is 
evident.  Solid  cloth  parachutes  inflate  more  rapidly  at  higher  altitudes  due  to  the 
decrease  in  air  density.  At  constant  trajectory  velocity  or  constant  dynamic  pressure 
the  effects  of  the  particular  value  of  n  converge  and  the  inflation  time  is  essentially 
independent  of  n  and  k.  Therefore,  Equation  (11)  for  high  altitude  horizontal 
trajectories  is  appropriate.  The  inflation  distance  versus  altitude  of  Figure  13  for 
n  =  0.500  is  the  same  whether  constant  velocity  or  dynamic  pressure  profile  is  used. 
The  realistic  n  =  0.632  of  Figure  14  gives  a  slight  variation  in  inflation  distance. 


APPLICATION  OF  CLOTH  PERMEABILITY  TO  THE  CALCULATION  OF  THE 
INFLATION  TIME  OF  PARACHUTES 

The  mass  outflow  through  the  pressurized  region  of  an  inflating  parachute  at 
any  instant  is  dependent  upon  the  instantaneous  canopy  area  which  is  subjected  to 
airflow  and  the  rate  of  airflow  through  that  area.  The  variation  of  pressurized  canopy 
surface  area  as  a  function  of  reference  time  ratio  t/to,  was  earlier  defined  as  being 
proportional  to  the  instantaneous  drag  area  ratio,  leaving  the  rate-of-airflow  problem 
to  solve.  The  permeability  parameter  of  cloth  was  a  natural  choice  for  determining 
the  rate  of  airflow  through  the  cloth  as  a  function  of  pressure  differential  across  the 
cloth.  Heretofore,  these  data  have  been  more  of  a  qualitative,  rather  than 
quantitative,  value.  A  method  of  analysis  was  developed  wherein  a  generalized  curve 
of  the  form  P  =  k(  AP)n  was  fitted  to  the  cloth  permeability  data  for  a  number  of 
different  cloths  and  gives  surprisingly  good  agreement  over  the  pressure  differential 
range  of  available  data.  The  pressure  differential  was  related  to  the  trajectory 
velocity  and  altitude  conditions  to  give  a  generalized  expression  which  can 
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CONSTANT  VELOCITY  PROFILE. 
POLAR  PLOT  OF  THE  RESULTS 
OF  EXAMPLE  1  OF  REFERENCE  5. 


ALTITUDE 

y 

(DEG) 

n  =  0.500 

n  =  0.662 

*0 

(SEC) 

*0 

(SEC) 

-90 

1.054 

1.416 

-75 

1.047 

1.401 

-60 

1.026 

1.362 

-45 

0.998 

1.31 1 

-30 

0.970 

1.259 

LU 

> 

-15 

0.940 

1.213 

LU 

—I 

0 

0.913 

1.171 

< 

LU 

15 

0.894 

1.138 

C/> 

30 

0.874 

1.108 

45 

0.860 

1.088 

60 

0.849 

1.074 

75 

0.841 

1.065 

90 

0.839 

1.062 

-90 

0.904 

1.171 

-  75 

0.900 

1.163 

-60 

0.894 

1.142 

-45 

0.876 

1.11 1 

-30 

0.854 

1.079 

LL 

-15 

0.833 

1.046 

O 

o 

0 

0.813 

1.018 

o 

CO 

15 

0.794 

0.992 

30 

0.779 

0.972 

45 

0.767 

0.955 

60 

0.758 

0.944 

75 

0.753 

0.938 

90 

0.751 

0.936 

-90 

0.513 

0.590 

-75 

0.512 

0.589 

-60 

0.509 

0.585 

-45 

0.504 

0.579 

-30 

0.499 

0.572 

h- 

LL 

-15 

0.493 

0.565 

o 

o 

0 

0.486 

0.559 

o 

o 

15 

0.481 

0.551 

CM 

30 

0.476 

0.545 

45 

0.471 

0.539 

60 

0.469 

0.537 

75 

0.467 

0.534 

90 

0  466 

0.534 

FIGURE  10.  EFFECTS  OF  TRAJECTORY  LAUNCH  ANGLE,  ALTITUDE,  AND  CLOTH  RATE  OF 
AIRFLOW  ON  THE  PARACHUTE  INFLATION  REFERENCE  TIME  AT  CONSTANT 
VELOCITY,  EXAMPLE  1  OF  REFERENCE 
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CONSTANT  DYNAMIC  PRESSURE 
PROFILE.  POLAR  PLOT  OF  THE 
RESULTS  OF  EXAMPLE  1  OF 


ALTITUDE 

FT 

SYMBOL 

n  =  0.5 

n  =  0.632 

0 

o 

• 

3,000 

□ 

■ 

20,000 

A 

▲ 

LD 

n  =  0.500 

n  =  0.662 

ALTITUD 

y 

(DEG) 

to 

(SEC) 

to 

(SEC) 

-90 

1.054 

1.418 

-75 

1.048 

1.400 

-60 

1.027 

1.363 

-45 

0.999 

1.312 

_ i 

-30 

0.968 

1.260 

LLi 

> 

-  15 

0.941 

1.214 

LLi 

—1 

0 

0.913 

1.172 

< 

LU 

15 

0.894 

1.138 

CO 

30 

0.874 

1  109 

45 

0.860 

1.089 

60 

0.850 

1.074 

75 

0.842 

1.066 

90 

0.840 

1.063 

-90 

0.862 

1  115 

-75 

0.859 

1.109 

-60 

0.849 

1.090 

-45 

0.834 

1  065 

-30 

0.814 

1  036 

h- 

LL 

-  15 

0.795 

1.008 

O 

o 

0 

0.778 

0.982 

o 

prs 

15 

0.762 

0.958 

30 

0.749 

0.941 

45 

0.738 

0.925 

60 

0.730 

0  916 

75 

0.725 

0  909 

90 

0  724 

0.908 

90 

0.366 

0  433 

75 

0.365 

0  433 

-60 

0.364 

0  432 

-45 

0.362 

0.429 

-30 

0  361 

0.428 

H 

LL 

-  15 

0,358 

0.425 

o 

o 

0 

0.357 

0.422 

o 

6 

15 

0.354 

0.419 

eg 

30 

0.352 

0.416 

45 

0.350 

0.414 

60 

0.349 

0  412 

75 

0.348 

0.41  1 

90 

0.348 

0.41  1 

FIGURE  11.  EFFECTS  OF  TRAJECTORY  LAUNCH  ANGLE,  ALTITUDE,  AND  CLOTH  RATE  OF 
AIRFLOW  ON  THE  PARACHUTE  INFLATION  REFERENCE  TIME  AT  CONSTANT 
DYNAMIC  PRESSURE,  EXAMPLE  1  OF  REFERENCE  5 
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0  0.2  0.4  0.6  0.8  1.0  1.2 


INFLATION  REFERENCE  TIME  (SEC) 


FIGURE  12.  EFFECT  OF  ALTITUDE  ON  THE  INFLATION  REFERENCE  TIME  "To" 
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0  100  200  300  400 


INFLATION  DISTANCE  (FEET) 

FIGURE  13.  EFFECT  OF  ALTITUDE  ON  THE  INFLATION  DISTANCE 
n  =  0.500,7  =  0\r  =  0 
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0  100  200  300  400 


INFLATION  DISTANCE  (FEET) 


FIGURE  14.  EFFECT  OF  ALTITUDE  ON  THE  INFLATION  DISTANCE 
N  =  0.632,  7  =  O',  t  =  0 
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be  used  in  the  finite,  intermediate,  and  infinite  Ballistic  Mass  Ratio  ranges.  This 
approach  varies  the  rate  of  airflow  for  the  particular  trajectory  conditions  at  each 
instant  during  inflation. 

Measured  and  calculated  permeability  pressure  data  for  several  standard  cloths 
are  illustrated  in  Figure  15.  This  method  has  been  applied  to  various  types  of  cloth 
between  the  extremes  of  a  highly  permeable  3-momme  silk  to  a  relatively  impervious 
parachute  pack  container  cloth  with  reasonably  good  results,  see  Figure  16. 

The  canopy  average  pressure  coefficient,  Cp,  is  defined  as  the  ratio  of  the 
pressure  differential  across  the  cloth  to  the  dynamic  pressure  of  the  free  stream. 


^  AP  P(internal)- P(external) 

p=  vT^2 


(25) 


where  V  is  the  instantaneous  trajectory  velocity. 

The  permeability  expression,  P  =  k(AP)n  becomes 


The  inflation  time  of  solid  cloth  parachutes  decreases  as  the  operational 
altitude  increases.  This  effect  can  be  explained  by  considering  the  ratio  of  the  mass 
outflow  through  a  unit  cloth  area  to  the  mass  inflow  through  a  unit  mouth  area. 

M'  =  mass  flow  ratio  =  mass  outflow 

mass  inflow 


where 

mass  outflow  =  P  p  slugs  (per  ft2  cloth  area) 

ft^  sec 


and 


mass  inflow  =  Vp  slugs  (per  ft2  inflow  area) 
ft12  sec 


Therefore,  the  mass  flow  ratio  becomes 

•r.-EE- 

p  V  V 


M' 


(26) 


Effective  porosity  C,  is  defined  as  the  ratio  of  the  velocity  through  the  cloth,  u, 
to  a  fictitious  theoretical  velocity,  v,  which  will  produce  the  particular  AP  =  pv2/2. 

effective  porosity,  C  =  u/v  (27) 
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FIGURE  15.  NOMINAL  POROSITY  OF  PARACHUTE  MATERIAL  VERSUS  DIFFERENTIAL 
PRESSURE 
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Comparison  of  the  mass  flow  ratio  and  previously  published  effective  porosity 
data  is  shown  in  Figure  17.  The  effects  of  altitude  and  velocity  on  the  mass  flow  ratio 
are  presented  in  Figures  18  and  19  for  constant  velocity  and  constant  altitude.  The 
decrease  of  cloth  permeability  with  altitude  is  evident.  The  reduced  outflow  rate 
causes  the  canopy  to  inflate  more  quickly  as  the  altitude  is  raised. 

The  relationship  of  k,  n,  and  altitude  on  the  inflation  reference  time  is 
illustrated  in  Figure  20. 

The  permeability  constants  "k"  and  "n"  can  be  determined  from  the 
permeability-pressure  differential  data  as  obtained  from  an  instrument  such  as  a 
Frazier  Permeameter.  Two  data  points,  "A"  and  "B,"  are  selected  in  such  a  manner 
that  point  "A"  is  in  a  low-pressure  zone  below  the  knee  of  the  curve,  and  point  "B"  is 
located  in  the  upper  end  of  the  high-pressure  zone,  as  shown  in  Figure  21. 

Two  measurements  at  1/2  inch  of  water  and  20  inches  of  water  appear  to  be  good 
data  points  if  both  are  available  on  the  same  sample.  Substituting  the  data  from 
points  "A"  and  "B"  into  P  =  k  (AP)n: 


UPa/ 


(28) 


k  = 


Pb 

(APB]n 


(29) 


Another  method  of  determining  k  and  n  is  to  adapt  a  least  squares  fit  through 
all  of  the  cloth  rate  of  airflow  -  pressure  differential  measurements.  A  true  average 
value  of  k  and  n  in  any  given  parachute  requires  a  large  number  of  permeameter  rate 
of  airflow  tests  due  to  the  variation  of  the  woven  cloth  airflow  rates.  Their  values 
may  also  vary  for  different  parachutes  in  the  same  lot.  This  is  most  likely  one  of  the 
causes  of  performance  variations  for  similar  test  conditions. 


PARACHUTE  AVERAGE  PRESSURE  COEFFICIENT 

The  average  pressure  coefficient,  Cpav,  of  a  parachute  is  the  parachute  drag 
coefficient  referenced  to  the  projected  canopy  area. 

The  average  pressure  coefficient  for  the  steady-state  parachute  was  determined 
as  follows: 
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FIGURE  17.  THE  EFFECTIVE  POROSITY  OF  PARACHUTE  MATERIALS  VERSUS 
DIFFERENTIAL  PRESSURE 
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FIGURE  18.  EFFECT  OF  ALTITUDE  ON  MASS  FLOW  RATIO  AT  CONSTANT  VELOCITY 


FIGURE  19.  EFFECT  OF  VELOCITY  ON  MASS  FLOW  RATIO  AT  CONSTANT  DENSITY 
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20  KFT 


INFLATION  REFERENCE  TIME  (SEC)  INFLATION  REFERENCE  TIME  (SEC) 


50  KFT  100  KFT 


INFLATION  REFERENCE  TIME  (SEC)  INFLATION  REFERENCE  TIME  (SEC) 


CLOTH 

k 

n 

MIL-C-7020,  TYPE  III 

1.46042 

0.63246 

MIL-C-8021,  TYPE  1 

0.99469 

0.57403 

MIL-C-802 1 ,  TYPE  II 

0.44117 

0.60492 

FIGURE  20.  EFFECT  OF  CANOPY  CLOTH  CONSTANTS  "k"  AND  "n”  ON  THE  PARACHUTE 
INFLATION  REFERENCE  TIME  "T0” 
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FIGURE  21.  LOCATION  OF  DATA  POINTS  FOR  DETERMINATION  OF  "k”  AND  "n” 


The  steady-state  drag  force,  in  a  wind  tunnel  for  example,  may  be  written 


f  -  qcDs0 

but  it  may  also  be  written 


F  = 


A  P  S 

av  p 


where  Sp  is  the  projected  area  of  the  inflated  canopy  and  APav  is  the 
average  pressure  differential  acting  on  Sp. 

^FavSp  F  =  q  CDS0 


(30) 


(30a) 


S  =  -7-  D2 
°  4  o 
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ratios  of  Dp/D0  from  Reference  6  are  presented  in  Table  2  for  several  types 
of  parachutes. 


TABLE  2.  RANGE  OF  AVERAGE  STEADY-STATE  CANOPY  PRESSURE  COEFFICIENTS 


PARACHUTE 

TYPE 

CD 

RANGE 

Dp/D0 

RANGE 

Cpav 

RANGE 

FLAT 

0.75 

0.70 

1.531 

CIRCULAR 

0.80 

0.67 

1.782 

EXTENDED 

0.78 

0.70 

1.592 

SKIRT 

0.87 

0.66 

1.997 

CROSS 

0.60 

0.72 

1.157 

0.78 

0.66 

1.791 

For  a  flat  circular  parachute  with  a  Cd  =  0.75  and  Dp/D0  =  0.67 

c  =  075 

Pav  (0.67)2 
Cn  =  1.671 

Mav 

The  inflation  geometry  as  a  function  of  t/to  is  constant  with  altitude.  It  can  be 
deduced  from  this  that  if  the  inflating  geometry  is  independent  of  altitude,  the  force 
distribution  which  causes  the  geometry  to  develop  must  also  be  independent  of 
altitude.  These  forces  come  from  the  pressure  distribution  along  the  gore  panel. 
Therefore,  at  any  given  t/to  &  definite,  repeatable  pressure  distribution  exists  which 
varies  with  t/to,  but  is  constant  for  all  altitudes  and  the  Cpav  is  constant  with  altitude. 
Figure  22  gives  the  effect  of  varying  the  magnitude  of  the  average  pressure 
coefficient  at  altitudes  from  sea  level  to  100,000  feet.  At  sea  level  the  magnitude  of 
the  average  pressure  coefficient  affects  the  value  of  the  inflation  reference  time.  The 
effect  is  magnified  as  the  exponent  V  increases.  As  the  altitude  increases  the 
particular  magnitude  of  the  average  pressure  coefficient  is  less  significant. 

Table  2  presented  the  range  of  Cpav  that  can  be  expected  for  a  definite 
parachute  type.  Figure  22  presents  the  effects  of  the  variation  of  the  average 
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AVERAGE  CANOPY  PRESSURE  COEFFICIENT 
DURING  INFLATION  INCLUDING  THE  VENT 


"n"  «  0.632 
Eq. (10) 


20  K  -  F \ 


50  K  -  FT 


100  K-  FT 


H 


c 


pav 


FIGURE  22.  EFFECT  OF  PRESSURE  COEFFICIENT  AND  ALTITUDE  ON  THE 
UNFOLDING  TIME 
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f*ressure  coefficient  on  the  inflation  reference  time  for  a  given  parachute  system.  The 
argest  inflation  reference  time  variation  occurs  at  sea  level.  As  the  altitude  is  raised 
the  variation  decreases  to  a  constant  time  value.  As  n  approaches  1  the  noted  effects 
are  accentuated. 

The  Air  Force  Flight  Dynamics  Laboratory  conducted  field  tests  of  a  28-foot  (D0) 
Solid  Flat  Circular  Parachute  (SFCP)  system  at  altitudes  of 6,000, 13,000,  and  21,000 
feet.  A  result  of  these  tests  was  the  demonstration  of  the  repeatability  of  the  canopy 
projected  area  to  surface  area  ratio  at  all  test  altitudes.  With  the  data  of  Figure  1,  the 
average  pressure  coefficients  for  the  28-foot  SFCP  can  be  estimated  during  the 
inflation  process. 


(31) 


for  0  <.  t/to  <.  0.3 


C 


0.0024  +  0.07843 


(32) 


for  0.3  <.  tAo  <.1.0 


C 


0.0117  +  0.0034  (118.4)^0 


(33) 


when  17  =  t  =  0 
for  0  <.  t/to  <_  0.3 


C 


(34) 


av  0.0024  +  0.07843 
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for  0.3  <.t/to  <.1.0 


0.0117  +  0.0034  (118.4)  t/lo 


(35) 


Figure  23  shows  the  estimated  variation  of  Cp  during  inflation  for  the  following 
conditions:  Cd  =  0.75;  n  =  0.0024;  to  =  tr,  and  the  two  variations  of  Sp/S0  as  a 
function  of  t/tf  in  Figure  1. 
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FIGURE  23.  ESTIMATED  STEADY-STATE  AVERAGE  PRESSURE  COEFFICIENT  FOR  THE 
INFLATING  28-FOOT  (D0)  SOLID  FLAT  CIRCULAR  PARACHUTE  OF  FIGURE  1 
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When  calculating  the  inflation  reference  time  of  solid  cloth  parachutes,  I  have 
achieved  good  results  using  a  constant  value  of  steady-state  average  pressure 
coefficient  in  the  inflation  time  equations. 
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CONCLUSIONS 


1.  An  inflating  parachute's  transient  geometric  shapes,  volumes,  and  drag  areas 
are  independent  of  altitude,  velocity,  and  system  mass  restraints. 

2.  The  parachute  inflation  distance  (Vsto)  axiom  has  been  verified  and 
demonstrated  to  be  a  constant  which  is  a  function  of  altitude. 

3.  The  parachute  inflation  time-opening  shock  force  axiom  has  been  demonstrated 
to  appear  to  be  true,  but  really  is  not. 

4.  The  opening  shock  force  of  a  given  parachute  system  is  determined  by  the 
initial  design  requirements,  type  of  parachute,  and  canopy  airflow,  as  well  as  the 
operational  altitude,  velocity,  and  trajectory  deployment  angle. 

5.  The  inflation  reference  times  of  solid  cloth  parachutes  vary  as  the  weight-to- 
drag-area  ratio  increases. 

6.  The  inflation  distance  formula  (12)  shows  that  the  performance  depends  on 
system  area  ratios,  mass  ratios,  canopy  airflow,  and  the  Ballistic  Mass  Ratio  scale 
parameter  per  foot  of  length. 

7.  The  inflation  distance  formula  effective  mouth  area  term  yields  some  clues 
about  critical  velocity. 

8.  Inflation  time  formulas  which  require  the  line  stretch  velocity  to  be  raised  to 
the  0.9  power,  can  achieve  the  same  result  by  dividing  the  line  stretch  velocity  by  2 
and  setting  the  exponent  of  Vs  to  unity. 

9.  Equation  (11)  shows  good  agreement  with  the  empirical  inflation  time  in  the 
low  altitude  regime,  and  shows  a  reduction  in  inflation  time  as  the  altitude  rises. 

10.  The  reduction  of  inflation  reference  time  and  distance  is  due  to  the  lowering  of 
air  density  as  the  altitude  increases. 

11.  The  cloth  rate  of  airflow  formula  is  suitable  to  high  rates  of  airflow  and  low 
rates  of  airflow  cloths. 

12.  The  mass  flow  ratio,  M’,  is  similar  to  effective  porosity,  C. 

13.  At  low  altitudes  the  particular  value  of  the  cloth  airflow  constants  have  a 
considerable  effect  on  the  inflation  reference  time.  As  the  altitude  increases,  the 
effect  of  k  and  n  on  the  inflation  reference  time  is  reduced.  At  an  altitude  of  100,000 
feet,  the  time  variation  is  minimal. 

14.  The  average  pressure  coefficient  of  a  parachute  is  the  parachute  drag 
coefficient  based  upon  the  projected  area  of  the  inflated  canopy.  The  intermediate 
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instantaneous  average  pressure  coefficients  can  be  determined  if  the  instantaneous 
drag  force  and  projected  area  are  known. 

15.  The  effect  of  the  average  pressure  coefficient  on  the  variation  of  the  inflation 
reference  time  is  most  pronounced  at  sea  level  and  least  pronounced  at  high  altitude. 
As  the  exponent  n  increases  from  0.5,  the  effects  are  intensified. 
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A  TECHNIQUE  FOR  THE  CALCULATION  C?  THE  OPENING-SHOCK  FORCES  FOR 
SEVERAL  TYPES  OF  SOLID  CLOTH  PARACHUTES 

W.  P.  Ludtke 

Laval  Ordr.aree  Laboratory 
Silver  Spring,  Maryland 


Abstract 


An  analytical  method  of  calculating 
parachute  opening-shock  forces  based  upon 
wind-tunnel  derived  drag  area  time  signa¬ 
tures  of  several  solid  cloth  parachute 
types  ir,  conjunction  with  a  seale  faetor 
and  retardation  system  steady- state  param¬ 
eters  has  been  developed.  Methods  of 
analyzing  the  inflation  time,  geometry, 
eloth  airflow  properties  and  materials 
elasticity  are  included.  The  effects  of 
mass  ratio  and  altitude  on  the  magnitude 
and  time  of  occurrence  of  the  maximum 
opening  shoek  are  consistent  with  observed 
field  test  phenomena. 

I.  Introduction 


In  1965 y  the  Haval  Ordnance  Laboratory 
(NOL)  was  engaged  in  a  project  which 
utilized  a  35-foot-diameter,  10-pereent 
extended-skirt  parachute  (type  T-10)  as 
the  second  stage  of  a  retardation  system 
for  a  250-pound  payload.  Deployment  of 
the  T-10  paraehute  was  to  be  accomplished 
at  an  altitude  of  100,000  feet.  In  this 
rarefied  atmosphere,  the  problem  was  to 
determine  the  seeor.d  stage  deployment 
conditions  for  successful  operation.  A 
seareh  of  available  field  test  ir formation 
indicated  a  lack  of  data  on  the  use  of 
solid  eloth  parachutes  at  altitudes  above 
30,000  feet. 

The  approach  to  this  problem  was  as 
follows :  Utilizing  existing  wind- tunnel 

data,  low-altitude  field  test  data,  and 
reasonable  assumptions,  a  unique  engi¬ 
neering  approach  to  the  inflation  time 
and  opening-shock  problem  was  evolved  that 
provided  satisfactory  results.  Basically, 
the  method  combines  a  wind-tunnel  derived 
drag  area  ratio  signature  as  a  function 
of  deployment  time  with  a  seale  faetor  and 
Newton*  s  second  lav;  of  motion  to  analyze 
the  velocity  and  foree  profiles  during 
deployment.  The  parachute  deployment 
sequence  is  divided  into  two  phases.  The 
first  phase,  called  "unfolding  phase," 
where  the  canopy  is  undergoing  changes  in 
shape,  is  considered  to  be  inelastie  as 
the  parachute  inflates  initially  to  its 
steady-state  aerodynamic  size  for  the 
first  time.  At  this  point,  the  "elastic 
phase"  is  entered  where  it  is  considered 
that  the  elasticity  of  the  paraehute 
materials  enters  the  problem  and  resists 
the  applied  forces  until  the  canopy  has 
reached  full  inflation. 

The  developed  equations  are  in  agree¬ 
ment  with  the  observed  performance  of 
solid  cloth  parachutes  in  the  field,  such 
as  the  deerease  of  inflation  time  as 


altitude  Increases,  effects  of  altitude  on 
opening- shock  foree,  fir  lie  a*~d  infinite 
mass  operation,  and  inflat io'  distance. 

II.  Development  of  Velocity  Ratio  and 
Foree  Ratio  Equations  During 
the  Unfolding  Phase  of 
Paraehute  Deployment 


The  parachute  deployment  v:ould  take 
place  in  a  horizontal  attitude  in  accord¬ 
ance  with  Newton's  seeond  lav;  of  motion. 


ZF  ~  ma 


1  c  _  w 

-  p  V  CnS  =  - 


dV 

at 


It  was  recognized  that  other  factors,  such 
as  included  air  mass,  apparent  mass,  and 
their  derivatives,  also  contribute  forces 
aeti*  z  on  the  system.  Since  definition 
of  tiiese  parameters  was  difficult,  the 
analysis  was  conducted  in  the  simplified 
form  shown  above.  Comparison  o?  calculated 
results  ar.d  test  results  indicated  that 
the  omitted  terms  have  a  small  effect. 


r  cD-  -  v  £ a 
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Multiplying  the  right-hand  side  of  cqua- 


Multiplylr.g 
tion  (1)  by 


y 
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and  rearranging 
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O  CDSo 
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In  order  to  integrate  the  left-hand  term 
of  equation  (2),  the  drag  area  ratio  must 
be  defined  for  the  type  of7 parachute  under 
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analysis  ac  a  function  of  deployment 

**•  .,f>p  fir*  r.  4 

LJ...  ,  '-q  • 


FIG  1  TYPICAL  INFINITE  MASS  FORCE  TIME  HISTORY  OF  A 
SOLID  CLOTH  PARACHUTE  IN  A  WIND  TUNNEL 
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REPRODUCEO  FROM  REFERENCE  (1) 


FIG.  3  TYPICAL  FORCE  TIME  CURVE  FOR  A  10%  EXTENDED 

SKIRT  PARACHUTE  UNOER  INFINITE  MASS  CONOITIONS. 


I 

150 « 


FIG.  4  TYPICAL  FORCE  TIME  CURVE  FOR  A  PERSONNEL 
GUIOE  SURFACE  PARACHUTE  UNOER  INFINITE 
MASS  CONOITIONS 


figure  1  illustrates  a  typical  solid 
clcth  parachute  wir.d -tunnel  infinite  mass 
force- time  history  after  snatch.  In 
infinite  mass  deployment ,  the  maximum  size 
arc  maximum  shock  force  occur  at  the  time 
of  full  inflation,  t^.  However,  t^  is 
inappropriate  for  analysis  since  it  is 
dependent  upon  the  applied  load,  structural 
strength,  and  materials  elasticity.  The 
reference  time,  t^,  where  the  parachute  has 
attained  its  steady-state  aerodynamic  size 
for  the  first  time,  is  used  as  the  basis  for 
performance  calculations. 

At  any  instant  during  the  unfolding 
phase,  the  force  ratio  F/Fs  can  be  deter¬ 
mined  as  a  function  of  the  time  ratio, 
t/-c  • 


7pVCDSo 


Since  the  wind-tunnel  velocity  and  density 
are  constant  during  infinite  mass  deployment 


FIG.  2  TYPICAL  FORCE  TIME  CURVE  FOR  A  SOLID  FLAT  PARA 
CHUTE  UNOER  INFINITE  MASS  CONOITIONS. 


FIG  5  TYPICAL  FORCE  TIME  SIGNATURE  FOR  THE  ELLIPTICAL 
PARACHUTE  UNDER  INFINITE  MASS  CONDITIONS 


TIME  (SEC) 


FIG  6  TYPICAL  FORCE-TIME  SIGNATURE  FOR  THE  RING 

SLOT  PARACHUTE  20%  GEOMETRIC  POROSITY  UNDER 
INFINITE  MASS  CONOITIONS 


Infinite  mass  opening- shock  signatures 
cf  several  types  of  parachutes  are  pre¬ 
sented  in  Figures  2  through  6.  Analysis 
of  nhese  signatures  usirg  the  force  ratio, 
F/Fc,  -  time  ratio,  t/tG,  technique 
indicated  a  similarity  in  the  performance 
of  the  various  sclid  cloth  types  of 
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parachutes  which  were  examined.  The  geo¬ 
metrically  porous  ring  slot  parachute 
displayed  a  completely  different  signature, 
as  was  expected.  These  data  are  illus¬ 
trated  in  Figure  7.  If  an  initial  boundary 
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FIG  7  DRAG  AREA  RATIO  VS  TIME  RATIO 


condition  of  CjjS/CjjSq  =  0  at  time  t/tQ  =  0 

is  assumed,  then,  the  data  can  be  approxi¬ 
mated  by  fitting  a  curve  of  the  form 


CDSo 


(3) 


At  the  tine  that  equation  (‘;)  was  ascer¬ 
tained,  it  suggested  that  the  geometry  of 
the  deploying  parachute  was  independent 
of  density  ana  velocity.  It  was  also 
postulated  that  although  this  expression 
had  been  determined  for  the  infinite  mass 
condition,  it  would  also  be  true  for 
the  finite  mass  case.  This  phenomenon 
has  since  been  independently  observed  and 
con firmed  oy  Berndt  and  2a  Tee se  in 
reference  (2). 

Since  the  drag  area  ratio  was  determined 
from  actual  parachute  deployments ,  it  was 
assumed  that  the  effects  of  apparent  mass 
and  included  mass  on  the  deployment  force 
history  were  accommodated. 

The  right-hand  term  of  equation  (2) 
contains  the  expression 

wv.Wo 


This  term  can  be  visualised  as  shown  in 
Figure  8  to  be  a  ratio  of  the  retarded 
mass  (including  the  parachute)  to  an 
associated  mass  of  atmosphere  contained 
in  a  right  circular  cylinder  which  is 
generated  by  moving  an  inflated  parachute 
of  area  0^fo  for  a  distance  equal  to  the 
product  of  Vst0  through  an  atmosphere  of 
density,  o. 


A  more  realistic  drag  area  ratio  expres¬ 
sion  was  determined  which  includes  the 
effect  of  initial  area  at  line  stretch. 


CDS 

CDSo 


-ft1-’)  (r)3  +  T  (t) 


FIG.  8  VISUALIZATION  OF  THE  MASS  RATIO  CONCEPT 


The  mass  ratio,  M,  is  the  scale  factor 
which  controls  the  velocity  and  force 
profiles  during  parachute  deployment. 
Substituting  K  ar.d  CpS/Cp$o  into  equa¬ 
tion  (2),  integrating,  and  solving  for 

V/V 
'  s 


where  r\  is  the  ratio  of  the  projected 
mouth  area  at  line  stretch  to  the  steady- 
state  projected  frontal  area.  Expanding 
equation  (4) 
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The  instantaneous  shock  factor  is  defined 
as 


x. 


L 


s 


s2cDSo 


Figures  9  and  10  illustrate  the  velocity 
and  force  profiles  generated  from  equations 
(7)  and  (8)  for  initial  projected  area 
ratios  of  r\  =  0,  and  0.2  with  various  mass 
ratios . 


IV.  Methods  for  Calculation  of 
the  Reference  Tine ,  tG 


If  the  altitude  variation  during  deployment 
is  small,  then,  the  density  may  be  consid¬ 
ered  as  constant 


from  equations  (5)  and  (7) 


III.  Maximum  Shock  Force  and  Time  of 
Occurrence  During  the  Unfolding  Phase 

The  time  of  occurrence  of  the  maximum 
instantaneous  shock  factor,  Xj_,  is  diffi¬ 
cult  to  determine  for  the  general  case. 
However,  for  ^  =  0,  the  maximum  shock 
factor  and  time  of  occurrence  are  readily 
calculated.  For  q  =  0 


The  ratio  concept  is  an  ideal  method  to 
analyze  the  effects  of  the  various  para-m¬ 
eters  on  the  velocity  and  force  profiles 
of  the  opening  parachutes;  however,  a 
means  of  calculating  tG  is  required  before 
specific  values  can  be  computed.  Methods 
for  computing  the  varying  mass  flow  into 
the  inflating  canopy  mouth,  the  varying 
mass  flow  out  through  the  varying  inflated 
canopy  surface  area,  and  the  volume  of 
air,  which  must  be  collected  during 

the  inflation  process  are  required. 

Figure  11  represents  a  solid  cloth- 
type  parachute  canopy  at  some  instant 
(8)  during  inflation.  At  any  given  instant, 
the  parachute  drag  area  is  proportional 
to  the  maximum  inflated  diameter.  Also, 
the  maximum  diameter  in  conjunction  with 
the  suspension  lines  determines  the  inflow 
mouth  area  (A-A)  and  the  pressurized 
canopy  area  (B-B-B).  This  observation 
provided  the  basis  for  the  following 
assumptions.  The  actual  canopy  shape  is 
of  minor  importance. 

a.  The  ratio  of  the  instantaneous 
mouth  inlet  area  to  the  steady-state 
mouth  area  is  in  the  same  ratio  as  the 
instantaneous  drag  area. 


x 


i 


Setting 
to  time 
t/t  at 


the  derivative  of  x^  with  respect 
equal  to  zero  and  solving  for 

xi  max 


(9) 


and  the  maximum  shock  factor  is 


(10) 


Tl  =  CDS 

CDSo 

b.  The  ratio  of  the  instantaneous 
pressurized  cloth  surface  area  to  the 
canopy  surface  area  is  in  the  same  ratio 
as  the  instantaneous  drag  area. 


c.  Since  the  suspension  lines  in  the 
ur.pressurized  area  of  the  canopy  are 
straight,  a  pressure  differential  has  not 
developed,  and,  therefore,  the  net  air¬ 
flow/  in  this  zone  is  zero. 

Based  on  the  foregoing  assumptions, 
the  mass  flow  equation  can  be  written 


Equations  (9)  and 
of  M  s  -1  (0.19), 

of  H,  the  maximum 
the  elastic  phase 


(10)  are  valid  for  value- 
sine  e  for  larger  values 

shock  force  occurs  in 
of  inflation. 


dm  =  m  inflow  -  m  outflow 


P  ~ 7  =  PVAu  -  pa sp 

dt 
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FIG.  9  EFFECT  OF  INITIAL  AREA  ANDMASS  RATIO  ON  THE 
SHOCK  FACTOR  AND  VELOCITY  RATIO  DURING  THE 
UNFOLDING  PHASE  FOR  *?  «  0. 


t/tQ  UNFOLDING  TIME  RATIO 


FIG.  10  EFFECT  OF  INITIAL  AREA  AND  MASS  RATIO  ON  THE 
SHOCK  FACTOR  AND  VELOCITY  RATIO  DURING  THE 
UNFOLDING  PHASE  FOR  ij-  0.2 
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FIG.  11  PARTIALLY  INFLATED  PARACHUTE  CANOPY 


From  equation  (3) 

CnS  /.  \6 

CT 


Vo  Uo/ 


;  for  t|  »  0 


From  equation  (7) 


i.ip) 

T«  Vo) 


;  n  =  o 


From  equation  ( 2o ) 


r2n 


Measured  values  of  n  indicate  a  data 
range  from  0.574  through  0.771*  A 
convenient  solution  to  the  reference  time 
equation  evolves  when  n  is  assigned  a 
value  of  1/2.  Integrating  equation  (13) 
and  using 


V  t  M  = 
s  o 


2W 


LET  K.  = 


gP-o 

2W 


CDSo 


Alo  "  ASok 


ar 


= 


i4w 


KPVsCDSo 


Kl  , 

e  -  1 


(14) 


Equation  (l4)  expresses  the  unfolding 
reference  time,  t0,  in  terms  of  mass, 
altitude,  snatch  velocity,  airflow  char¬ 
acteristics  of  the  cloth,  and  the  steady- 
state  parachute  geometry.  Note  that  the 
term  gpVQ/W  is  the  ratio  of  the  included 
air  mass  to  the  mass  of  the  retarded 
hardware.  Multiplying  both  sides  of 
equation  (14)  by  Vs  demonstrates  that 


V  t  =  a  constant  which  is 

a  function  of  altitude 


Figures  12  and  13  indicate  the  para¬ 
chute  unfolding  time  and  unfolding 
distance  for  values  of  n  =  1/2  and  n  = 
0.63246.  Note  the  variation  and  conver¬ 
gence  with  rising  altitude.  The  opening- 
shock  force  is  strongly  influenced  by 
the  inflation  time.  Because  of  this,  the 
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value  of  tQ  calculated  by  using  a  real¬ 
istic  value  of  n  should  be  used  in  the 
lower  atmosphere. 

As  an  example  of  this  method  of  opening- 
shock  analysis,  let  us  examine  the  effect 
of  altitude  on  the  opening-shock  force  of 
a  T-10-type  parachute  retarding  a  '^00- 
pound  weight  from  a  snatch  velocity  of 
Vs  a  400  feet  per  second  at  sea  level. 
Conditions  of  constant  velocity  ar.d  con¬ 
stant  dynamic  pressure  are  investigated. 

The  results  are  presented  in  Figure  14. 

At  low  altitudes,  the  opening-shock  force 
is  less  than  the  steady-state  drag  force: 
however,  as  altitude  rises,  the  opening 
shock  eventually  exceeds  the  steady-state 
drag  force  at  some  altitude.  This  trend 
is  in  agreement  with  field  test  observations. 

V.  Correction  of  tQ  for 
Initial  Area  Effects 


The  unfolding  reference  time,  tQ, 
calculated  by  the  previous  methods  assumes 
that  the  parachute  inflates  from  zero  drag 
area.  In  reality,  a  parachute  has  a  drag 


UNFOLDING  TIME  ISEC.) 


FIG.  12  EFFECT  OF  ALTITUDE  ON  THE  UNFOLDING  TIME 
~to"  AT  CONSTANT  VELOCITY  AND  CONSTANT 
DYNAMIC  PRESSURE  FOR  n  -  1/2  AND  n  -  0.63296 


CLOTH  MIL  C  7020  TYPE  III 

35  FOOT,  10°'.  EXTENDED  SKIRT  CANOPY 


FIG.  13  EFFECT  OF  ALTITUDE  ON  THE  UNFOLDING  DIS¬ 
TANCE  AT  CONSTANT  VELOCITY  AND  CONSTANT 
DYNAMIC  PRESSURE  FOR  n  -  1/2  AND  n  -  0  63246 


FIG.  14  VARIATION  OF  STEADY  STATE  DRAG,  F  AND  MAXI¬ 
MUM  OPENING  SHOCK  WITH  ALTITUDE  FOR  CONSTANT 
VELOCITY  AND  CONSTANT  DYNAMIC  PRESSURE 


area  at  the  beginring  of  inflation.  Once 
tQ  has  been  calculated,  a  correction  can 
be  applied,  based  upon  what  Is  known  about 
the  initial  conditions. 
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Case  A  -  When  the  initial  projected 

area  is  known 


than  the  maximum  extensibility,  e 
the  materials.  r 


of 


Z'&j 


■■■(?) 


i  \V3 


calculated 


Development  of  the  analysis  in  the 
elastic  phase  of  inflation  is  similar  to 
the  technique  used  in  the  unfolding  phase. 
Newt  or. 1  s  second  lav/  of  motion  is  used, 
together  with  the  drag  area  ratio  signa¬ 
ture  and  mass  ratio 


=  /_l  y 

CDSo  VJ 


correc 


■f'-pr 

ted  L  \AC  / 


calculated 


(15) 


Caae  B  -  When  the  initial  drag  area 

ts  known 


fsfi  =  fl±\( 

CDSo  ^o/ 


/C  S  \1/6 

(_ ¥±\ 

\CDSo/ 


calculated 


which  is  still  valid,  as  shown  in  Figure  7 


=U,fe)‘ v 

to 


Integrating  and  solving  for  — 


where  —  is  the  velocity  ratio  of  the 


(17) 


corrected 


#  C  S  \ 

,(^l\ 

Wo  f 


The  mass  ratio  should  ^ow  be  adjusted  for 
the  corrected  tG  before  velocity  and  force 
profiles  are  determined. 

VI.  Opening-Shock  Force,  Velocity  Ratio, 
and  Inflation  Time  During  the 
Elastic  Phase  of  Parachute  Inflation 

The  mass  ratio,  M,  is  an  important 
parameter  in  parachute  analysis.  For 
values  of  M  <<  4/21,  the  maximum  opening- 
shock  force  occurs  early  in  the  inflation 
process,  and  the  elastic  properties  of 
the  canopy  are  not  significant.  As  the 
mass  ratio  approaches  <A  =  4/21,  the 
magnitude  of  the  opening-shock  force 
increases,  and  the  time  of  occurrence 
happens  later  in  the  deployment  sequence. 
For  mass  ratios  M  >  4/21,  the  maximum 
shock  force  will  occur  after  the  reference 
time,  to-  Parachutes  designed  for  high 
mass  ratio  operation  must  provide  a 
structure  of  sufficient  constructed 
strength,  Fc,  so  that  the  actual  elon¬ 
gation  of  the  canopy  under  load  is  less 


unfolding  process  at  time  t  «»  tQ. 


"calculated 

(16) 


s  1  +  ± 


il  -  .i)g  +  nU  -  Ti)  +  „2 


(18) 


The  instantaneous  shock  factor  in  the 
elastic  phase  becomes 


x  =_dl/jlY 

1  CDSo\Vs/ 

fe)6 


+  -1- 
7M 


(r)- 


(19) 


The  end  point  of  the  inflation  process 
depends  upon  the  applied  ’loads,  elasticity 
of  the  canopy,  and  the  constructed  strength 
of  the  parachute.  A  linear  load  elongation 
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relationship  is  utilized  to  determine  the 
maximum  drag  area. 


F 


max 


The  next  force  in  the  series  at  constant  q 
F,  =  X  F  — 

1  °  S  r 

kc 


F  e 


max 


(20) 


v/nere 


r  =  u  +  *0) 


The  force,  F,  is  initially  the  instanta¬ 
neous  force  at  the  end  of  the  unfolding 
process 


=  X  F 
o  s 


(21) 


where  XQ  is  the  shock  factor  of  the 
unfolding  phase  at  t  =  tQ 


Xo  =- 


1  + 


i  rn^if + +  „?] 

M  L  7  2  J 


(22) 

Since  the  inflated  shape  is  defined,  the 
drag  coefficient  is  considered  to  be 
constant,  and  the  instantaneous  force  is 
proportional  to  the  dynamic  pressure  and 
projected  area.  The  maximum  projected 
area  would  be  developed  if  the  dynamic 
pressure  remained  constant  during  the 
elastic  phase.  Under  very  high  mass 
ratios,  this  is  nearly  the  case  over  this 
very  brief  time  period;  but  as  the  mass 
ratio  decreases,  the  velocity  decay  has 
a  more  significant  effect.  The  simplest 
approach  for  all  mass  ratios  is  to 
determine  the  maximum  drag  area  of  the 
canopy  as  if  elastic  inflation  had 
occurred  at  constant  dynamic  pressure. 

Then  utilizing  the  time  ratio  determined 
as  an  end  point,  intermediate  shock 
factors  can  be  calculated  from  equation 
(19)  and  maximum  force  assessed. 

The  initial  force,  X0FS,  causes  the 
canopy  to  increase  in  projected  area.  The 
new  projected  area  in  turn  increases  the 
total  force  on  the  canopy  which  produces 
a  secondary  projected  area  increase.  The 
resulting  series  of  events  are  resisted 
by  the  parachute  materials.  The  parachute 
must,  therefore,  be  constructed  of  suffi¬ 
cient  strength  to  prevent  the  elongation 
of  the  materials  from  exceeding  the 
maximum  elongation. 


X  F 
o  s 


max 


(23) 


Subsequent  elongations  ir  the  system  can 
be  shown  to  be 


:1  =  eo  +  eo >‘ 


2v2 


€2  =  eo  U  +  eo  t1  +  €o)  ) 

The  required  canopy  constructed  strength 
can  be  determined  for  a  given  set  of 
deployment  conditions.  The  limiting  value 
of  the  series  (g^)  determines  the  end 
point  time  ratio. 


/M6,  CDSir.ax  _ 
VV  ”•  CDSo 


=  (1  +  H)‘ 


/M  =  /VwW- 

\V  V  CDSo  ' 


=  (1  +  e,J 


1/3 


(24) 


Figure  15  illustrates  the  r.aximum  drag 
area  ratio  as  a  function  of  gq. 


F*G  15  MAXIMUM  DRAG  AREA  RATIO  VS.  INITIAL 
ELONGATION 
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VII.  Application  of  Cloth  Permeability 
to  the  Calculation  of  the  Inflation 
Time  of  Solid  Cloth  Parachutes 


The  mass  outflow  through  the  pressur¬ 
ized  region  of  an  inflating  solid  cloth 
parachute  at  any  instant  is  dependent  upon 
the  canopy  area  which  is  subjected  to  air¬ 
flow  and  the  rate  of  airflow  through  that 
area.  The  variation  of  pressurized  area  as 
a  function  of  reference  time,  tQ,  was 
earlier  assumed  to  be  proportional  to  the 
Instantaneous  drag  area  ratio,  leaving  the 
rate- of- airflow  problem  to  solve.  The 
permeability  parameter  of  cloth  was  a 
natural  choice  for  determining  the  rate  of 
airflow  through  the  cloth  as  a  function  of 
pressure  differential  across  the  cloth. 
Heretofore,  these  data  have  been  more  of 
a  qualitative,  rather  than  quantitative, 
value.  A  new  method  of  analysis  v/as 
developed  wherein  a  generalized  curye  of 
the  form  P  =  k(AP)n  was  fitted  to  cloth 
permeability  data  for  a  number  of  different 
cloths  and  gives  surprisingly  good  agree¬ 
ment  over  the  pressure  differential  range 
of  available  data.  The  pressure  differ¬ 
ential  was  then  related  to  the  trajectory 
conditions  to  give  a  generalized  expression 
which  can  be  used  in  the  finite  mass  ratio 
range,  as  'well  as  the  infinite  mass  case. 
The  permeability  properties  were  trans¬ 
formed  Into  a  mass  flow  ratio,  Mf,  which 
shows  agreement  with  the  effective  porosity 
concept. 

Measured  and  calculated  permeability 
pressure  data  for  several  standard  cloths 
are  illustrated  in  Figure  1 6.  This  method 
has  been  applied  to  various  types  of  cloth 
between  the  extremes  of  a  highly  permeable 
3-momme  silk  to  a  relatively  Impervious 
parachute  pack  container  cloth  with 
reasonably  good  results,  see  Figure  17. 

The  canopy  pressure  coefficient,  Cp, 
is  defined  as  the  ratio  of  the  pressure 
differential  across  the  cloth  to  the 
dynamic  pressure  of  the  free  stream. 


Ap  _  P(lnternal)  P(external) 
q  1/2  pV2 


(25) 

where  V  is  based  on  equation  (7). 

The  permeability  expression,  P  =  k(AP)n 
becomes 

v2  n 

?  =  k(Cp  (2^) 


Although  some  progress  has  been  made 
by  Melzig  and  others  on  the  measurement  of 
the  variation  of  the  pressure  coefficient 
on  an  actual  inflating  canopy,  this 
dimension  an d  its  variation  with  time  are 
still  dark  areas  at  the  time  of  this 
writing.  At  the  present  time,  a  constant 
average  value  of  pressure  coefficient  is 


used  in  these  calculations.  Figure  lB 
presents  the  effect  of  pressure  coefficient 
and  altitude  on  the  unfolding  time  for 
constant  deployment  conditions. 

It  is  well  known  that  the  inflation 
time  of  solid  cloth  parachutes  decreases 
as  the  operational  altitude  increases. 

This  effect  can  be  explained  by  consid¬ 
ering  the  ratio  of  the  mass  outflow 
through  a  unit  cloth  area  to  the  mass 
inflow  through  a  unit  mouth  area. 


M* 


mass  flow  ratio 


mass  outflow 
mass  inflow 


where 


mass  outflow 


P.  (per  n2  cloth  area) 

°  ft  -sec 


and 


mass  Inflow  =  V  -  (per  ft2  inflow  area) 

P  ft  -sec 


REPRODUCED  FROM  REFERENCE  {4) 


FIG.  16  NOMINAL  POROSITY  OF  PARACHUTE  MATERIAL 
VS  DIFFERENTIAL  PRESSURE. 
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Therefore,  the  mass  flow  ratio  becomes 


Effective  porosity,  C,  is  defined  as 
the  ratio  of  the  velocity  through  the 
cloth,  u,  to  a  fictitious  theoretical 
velocity,  v,  which  will  produce  the  par¬ 
ticular  A  P=  l/2ov2. 


AVERAGE  CANOPY  PRESSURE  COEFFICIENT 
DURING  INFLATION  INCLUDING  THE  VENT 


FIG.  18  EFFECTOF  PRESSURE  COEFFICIENT  AND  ALTITUDE 
ON  THE  UNFOLDING  TIME. 


effective  porosity,  C  *  —  (28) 


Comparison  of  the  mass  flow  ratio  and 
previously  published  effective  porosity 
data  is  shown  in  Figure  19.  The  effects 
of  altitude  and  velocity  on  the  mass  flow 
ratio  are  presented  in  Figures  20,  and  21 
for  constant  velocity  and  constant 
altitude.  The  decrease  of  cloth  perme¬ 
ability  with  altitude  is  evident. 

The  permeability  constants  "k"  and  "n" 
can  be  determined  from  the  permeability 
pressure  differential  data  as  obtained 
from  an  instrument  such  as  a  Frazier 
Permeameter.  Two  data  points,  rA"  and 


FIG.  17  COMPARISON  OF  MEASURED  AND  CALCULATED 

PERMEABILITY  FOR  RELATIVELY  PERMEABLE  AND 
IMPERMEABLE  CLOTHS 


REPRODUCED  FROM  REFERENCE  (4) 

FIG.  19  THE  EFFECTIVE  POROSITY  OF  PARACHUTE 
MATERIALS  VS.  DIFFERENTIAL  PRESSURE 
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"B,"  are  selected  ir.  such  a  manner  that 
point  "A"  is  in  a  lev-pressure  zone  below 
the  knee  of  the  cur  e,  and  point  ,fBn  is 
located  in  the  upper  end  of  the  high- 
press  uro  zone,  as  shown  in  Figure  22. 

The  two  standard  measurements  of  1/2 
inch  of  water  and  20  inches  of  water 
appear  to  be  good  data  points  if  both  are 


FIG.  20  EFFECT  OF  ALTITUDE  ON  MASS  FLOW  RATIO  AT 
CONSTANT  VELOCITY 


VELOCITY  (F  P.S.) 


FIG  21  EFFECTOF  VELOCITY  ON  MASS  FLOW  RATIO  AT 
CONSTANT  DENSITY 


FIG  22  LOCATION  OF  DATA  POINTS  FOR  DETERMINATION 
OF  "k"  AND  “n" 


available  on  the  same  sample.  Substi¬ 
tuting  the  data  from  points  "A”  and  ,fBn 


into  P  =  k(AP)n  : 

.  tn(^) 


(?9) 


(30) 


VIII.  Determination  of  the 
Parachute  Included  Volume 
and  Associated  Air  Mass 


Before  the  reference  time,  tG,  and 
inflation  time,  tf,  can  be  calculated, 
the  volume  of  atmosphere,  Vo,  which  is 
to  be  collected  during  the  inflation 
process  must  be  accurately  known.  This 
requirement  dictates  that  a  realistic 
inflated  canopy  shape  and  associated 
volume  of  atmosphere  be  determined. 

Figure  23  was  reproduced  from  reference 
(5).  The  technique  of  using  lampblack 
coated  plates  to  determine  the  airflow 
patterns  around  metal  models  of  inflated 
canopy  shapes  was  used  by  the  investigator 
of  reference  (5)  to  study  the  stability 
characteristics  of  contemporary  parachutes, 
i.e.,  1943.  A  by-product  of  this  study 
is  that  it  is  clearly  shown  that  the 
volume  of  air  within  the  canopy  bulges 
out  of  the  canopy  mouth  (indicated  by 
arrov:s)  and  extends  ahead  of  the  canopy 
hem.  This  volume  must  be  collected  during 
the  inflation  process.  Another  neglected, 
but  significant,  source  of  canopy  volume 
exists  in  the  billowed  portion  of  the  gore 
panels . 


A  13 


NSWC  TR  88-292 


HEMISPHERE 


VENT  PARACHUTE 

REPRODUCED  FROM  REFERENCE  (5) 

FIG.  23  AIRFLOW  PATTERNS  SHOWING  AIR  VOLUME 
AHEAD  OF  CANOPY  HEM 


The  steady- state  canopy  shape  has  been 
observed  in  wind-tunnel  and  field  tests  to 
be  elliptical  in  profile.  Studies  of  the 
inflated  shape  and  included  volume  of 
several  parachute  types  (flat  circular* 

10  percent  extended  skirt*  elliptical* 
hemispherical*  ring  slot*  ribbon*  and 
cross)  are  documented  in  references  (6) 
and  (7).  These  studies  demonstrated  that 
the  steady-state  profile  shape  of  inflated 
canopies  of  the  various  types  can  be 
approximated  to  be  two  ellipses  of  common 
major  diameter*  2a*  and  dissimilar  minor 
diameters*  b  and  b1*  as  shown  in  Figure  24. 
It  was  also  shown  that  the  volume  of  the 
ellipsoid  of  revolution  formed  by  revolving 
the  profile  shape  about  the  canopy  axis 
v;as  a  good  approximation  of  the  volume  of 
atmosphere  to  be  collected  during  canopy 
inflation  and  included  the  air  volume 
extended  ahead  of  the  parachute  skirt  hem 
together  with  the  billowed  gore  volume. 


— 1  BILLOWED  GORE 
v!  v  CONTOUR 


REPRODUCED  FROM  REFERENCE  (6) 

FIG.  24  PARACHUTE  CROSS  SECTION  NOMENCLATURE 

Tables  I  and  II  arc  summaries  of  test 
results  reproduced  from  references  (6)  and 
(7)*  respectively*  for  the  convenience  of 
the  reader. 
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TABLE  I  SUMMARY  OF  PARACHUTE  SHAPE  TEST  RESULTS 
FOR  12-GORE  AMD  16-GORE  CONFIGURATIONS 


Parachute 

Type 

No.  of 

Gores 

Suspension 
Line  Length 
inches 

Velocity 

mph  fps 

2  a 

°o 

Scale  Factor,  K 
23  21 

°F  °B 

2a 

L 

N 

a 

Axes  Ratio 

b  b' 

a  a 

b  b' 

a  i 

Volume  in 

VH  Vc 

3 

V 

o 

Ol  X 

>  l> 

Flat  Circular 

12 

34 

50 

73 

645 

.650 

856 

6115 

8817 

1.4932 

4476 

4481 

6980 

1  56 

16 

34 

50 

73 

663 

.669 

.820 

.5558 

.9039 

1.4597 

4450 

4100 

7325 

1  65 

10%  Extended 

12 

34 

100 

147 

.663 

652 

881 

6424 

8860 

1  5284 

3928 

4400 

6783 

1.73 

Skirt 

16 

34 

17 

25 

654 

640 

785 

5580 

8502 

1  4082 

4051 

3920 

6197 

1.53 

Elliptical 

12 

34 

75 

110 

916 

.812 

5626 

.9657 

1  5283 

3322 

5405 

16 

34 

17 

25 

.875 

.800 

6169 

8163 

1  4332 

2726 

4405 

Hemispherical 

12 

34 

125 

183 

996 

1  254 

1.0005 

9080 

1  9085 

6224 

8666 

16 

34 

75 

110 

994 

1  185 

9129 

.9380 

1  8509 

5921 

8370 

Ringslot 

12 

34 

25 

37 

607 

654 

853 

.6566 

.8735 

1.530 

3800 

3650 

5903 

1  55 

16%  Geometric 

12 

34 

100 

147 

.616 

663 

922 

.6566 

8735 

1  530 

3800 

4198 

6166 

1  62 

Porosity 

12 

34 

200 

293 

.637 

686 

918 

6566 

8735 

1  530 

3800 

4624 

6826 

1  90 

16 

34 

25 

37 

.611 

.658 

827 

6004 

8890 

1  4894 

3800 

3763 

5685 

1  50 

16 

34 

100 

147 

.617 

664 

.864 

6004 

.8890 

1.4894 

3800 

3985 

6030 

1  59 

16 

34 

200 

293 

.645 

695 

.844 

6004 

8890 

1  4894 

3800 

4430 

6897 

1  82 

Ribbon 

12 

34 

25 

37 

586 

632 

859 

6558 

.8768 

1.5326 

3800 

3323 

5335 

1  40 

24%  Geometric 

12 

34 

100 

147 

.615 

.663 

.837 

6558 

8768 

1.5326 

3800 

3714 

6163 

1  62 

Porosity 

12 

34 

200 

293 

.632 

.681 

877 

6558 

.8768 

1.5326 

3800 

4280 

6683 

1  76 

16 

34 

25 

37 

.603 

650 

797 

5570 

8578 

1  4148 

3800 

3438 

5358 

1  41 

16 

34 

100 

147 

626 

.674 

791 

.5570 

8578 

1.4148 

3800 

3804 

5983 

1  57 

16 

34 

200 

293 

.648 

.698 

.781 

5570 

8578 

1.4148 

3800 

4164 

6656 

1.75 

Cross  Chute 

34 

25 

37 

.710 

543 

1.242 

.8867 

1  2776 

2  1643 

1928 

3768 

5798 

3  01 

w/L  -  264 

34 

100 

147 

707 

.540 

1  270 

.8867 

1.2776 

2  1643 

1928 

3810 

5712 

2  96 

34 

200 

293 

.716 

.547 

1.285 

8867 

1  2776 

2  1643 

1928 

4212 

5925 

307 

47 

25 

37 

.759 

580 

1.113 

.8494 

1  2512 

2  1006 

1928 

4052 

6868 

3  56 

47 

100 

147 

729 

557 

1  205 

8494 

1  2512 

2.1006 

1928 

3973 

5958 

3  09 

47 

200 

293 

775 

592 

1.110 

8494 

1  2512 

2.1006 

1928 

4292 

7303 

3  79 

REPRODUCED  FROM  REFERENCE  (6) 


TABLE  II  SUMMARY  OF  PARACHUTE  SHAPE  TEST  RESULTS 
FOR  24-GORE  AND  30-GORE  CONFIGURATIONS 


Parachute 

Type 

No  of 

Gores 

Suspension 
Line  Length 
inches 

Velocity 
mph  fps 

Scale  Factor  K 

2a  25 

O  Dc 

o  F 

N 

a 

b 

a 

Axes  Ratio 

b'  b  b 

a  a  +  i 

Volume  in  3 

V  V 

H  c 

vo 

V 

O 

Flat  Circulars 

24 

34 

50 

73 

.677 

679 

795 

.5758 

.8126 

1.3884 

4362 

4695 

7273 

1.67 

30 

34 

17 

25 

668 

669 

.827 

6214 

7806 

1  4020 

4342 

4626 

7027 

1.62 

10%  Extended* 

24 

34 

100 

147 

.665 

.648 

834 

5949 

.8771 

1.4720 

4138 

4446 

6930 

1  67 

Skirt 

30 

34 

17 

25 

650 

633 

825 

6255 

7962 

1.4127 

4172 

4076 

6265 

1.50 

Ring  Slot  16% 

24 

34 

25 

37 

663 

665 

824 

.5800 

.9053 

1.4853 

3591 

3878 

6031 

1  68 

Geometrically 

24 

34 

100 

147 

.680 

682 

819 

5800 

.9053 

1  4853 

3591 

4079 

6510 

1  81 

Porous 

24 

34 

200 

293 

.694 

696 

809 

5800 

.9053 

1  4853 

3591 

4270 

6924 

1  93 

30 

34 

25 

37 

677 

678 

788 

5800 

.9053 

1.4853 

3582 

3826 

6404 

1  79 

30 

34 

100 

147 

.684 

685 

802 

5800 

.9053 

1  4853 

3582 

4023 

6588 

1  84 

30 

34 

200 

293 

698 

699 

800 

.5800 

.9053 

1  4853 

3582 

4260 

7012 

1.96 

Ribbon  24% 

24 

34 

25 

37 

.671 

.673 

770 

5980 

.8187 

1  4167 

3591 

3591 

5968 

1  66 

Geometrically 

24 

34 

100 

147 

676 

678 

.813 

.5980 

8187 

1.4167 

3591 

3927 

6097 

1  70 

Porous 

24 

34 

200 

293 

687 

689 

.804 

5980 

.8187 

1  4167 

3591 

4061 

6389 

1  78 

30 

34 

25 

37 

655 

657 

.782 

.6021 

.8463 

1  4484 

3582 

3396 

5666 

1  58 

30 

34 

100 

147 

.669 

670 

784 

6021 

.8463 

1  4484 

3582 

3622 

6022 

1.68 

30 

34 

200 

293 

677 

679 

.823 

6021 

8463 

1  4484 

3582 

4002 

6256 

1  75 

'Since  this  parachute  was  "breathing  '  during  the  test,  several  photographs  were  taken  at  each  speed.  The  data  were  reduced  from  the  photograph  which  most 
reasonably  appeared  to  represent  the  equilibruim  state 


REPRODUCED  FROM  REFERENCE  (7) 
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A 


X.  List  of  Symbols 


t 


Steady-state  projected  area  of  the 
Inflated  parachute,  ft2 


-  Cloth  permeability  -  rate  of  air¬ 
flow  through  a  cloth  at  an  arbi¬ 
trary  differential  pressure, 

ft3/ft2/sec 


-  Instantaneous  canopy  mouth  area,  ft2  q 


-  Dynamic  pressure,  \t/ft2 


'Vto 


a 

2a 


-  Steady-state  inflated  mouth  area,  S  “A  ~  Instantaneous  inflated  canopy 

rt?  8  surface  area,  ft* 

-  Acceleration,  ft/sec2  s0°AS0  "  Canqpy  8Urface  area,  ft* 

-  Maximum  inflated  parachute  diameter 

of  gore  mains tarn,  ft  t  -  Instantaneous  time,  sec 


b  -  Minor  axis  of  the  ellipse  bounded 

by  the  major  axis  (ZS)  and  the  vent 
of  the  canopy,  ft 

b*  -  Minor  axis  of  the  ellipse  which 
includes  the  skirt  hem  of  the 
canopy,  ft 

C  -  Effective  porosity 

CD  -  Parachute  coefficient  of  drag 

C  -  Parachute  pressure  coefficient, 

p  relates  internal  and  external 

pressure  (aP)  on  canopy  surface  to 
the  dynamic  pressure  of  the  free 
3 1  re  am 

Dn  -  Nominal  diameter  of  the  aerodynamic 
decelerator  =  ft 

F  -  Instantaneous  force,  lbs 

F  -  Steady- state  drag  force  that  would 
be  produced  by  a  fully  epen  para¬ 
chute  at  velocity  Va,  lba 

F  -  Constructed  strength  of  the  para¬ 
chute,  lb3 

-  Maximum  opening- shock  force,  lbs 

max 

p 

g  -  Gravitational  acceleration,  ft/sec 

k  -  Permeability  constant  of  canopy 

cloth 


t  -  Reference  time  when  the  parachute 
has  reached  the  design  drag  area 
for  the  first  time,  sec 

t_  -  Canopy  inflation  time  when  the 
inflated  canopy  has  reached  its 
maximum  physical  site,  sec 

u  -  Air  velocity  through  cloth  in 
effective  porosity,  ft/sec 

v  -  Fictitious  theoretical  velocity 

used  in  effective  porosity,  ft/sec 

V  -  Instantaneous  system  velocity,  ft/eec 

V  “  System  velocity  at  the  time  t  -  t  , 

°  ft/sec  ° 

V  -  System  velocity  at  the  end  of 

suspension  line  stretch,  ft/eec 

V  -  Volume  of  air  which  must  be  collec¬ 

ted  during  the  inflation,  process, 

rt3 

W  -  Hardware  weight,  lb 

X£  -  Instantaneous  shock  factor 

X_  -  Shock  factor  at  the  time  t  -  t 

o  o 

p  -  Air  density,  slugs/ft^ 

T|  -  Ratio  of  parachute  projected  mouth 
area  at  line  stretch  to  the  steady- 
state  projected  area 


m  -  Mass,  slugs 

M  -  Mass  ratio  -  ratio  of  the  mass  of 
the  retarded  hardware  {including 
parachute)  to  a  mass  or  atmosphere 
contained  in  a  right  circular 
cylinder  of  length  (Vfit0),  face  area 
(CDS0),  and  density  (o) 

Mf  -  Maes  flow  ratio  -  latio  of  atmo¬ 
sphere  flowing  through  a  unit  cloth 
area  to  the  atmosphere  flowing 
through  a  unit  inlet  area  at 
arbitrary  pressure 


c  -  Instantaneous  elongation 

«  -  Maximum  elongation 

max 

e  -  Initial  elongation  at  the  beginning 

of  the  elastic  phase  of  inflation 

S.F.  -  Parachute  safety  factor  -  Fc/Fnax 


n  -  Permeability  constant  of  canopy 
cloth 
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Appendix  B 


A  GUIDE  FOR  THE  USE  OF  APPENDIX  A 


At  first  reading,  Appendix  A  may  appear  to  be  a  complicated  sys 
tern  of  analysis  because  of  the  many  formulae  presented.  Actually, 
once  understood,  the  technique  is  straightforward  and  uncomplicated 
The  author  has  attempted  to  simplify  the  algebra  wherever  possible. 
This  appendix  presents,  in  semi-outline  form,  a  guide  to  the 
sequence  of  calculations  because  the  analysis  does  require  use  of 
formulae  from  the  text,  not  necessarily  in  the  order  in  which  they 
were  presented.  Also,  the  user  can  be  referred  to  graphs  of  per¬ 
formance  to  illustrate  effects. 

In  order  to  compute  tQ,  other  parameters  must  be  obtained  from 
various  sources. 

I.  Determine  System  Parameters 

1.  CDS0,  drag  area,  ft^  obtained  from  design  requirement. 

2.  Vs,  fps,  velocity  of  system  at  suspension  line  stretch. 

3.  p,  slugs/ft-^,  air  density  at  deployment  altitude. 

4.  W,  lb,  system  weight  (including  weight  of  the  parachute) 
from  design  requirements. 

5.  VD,  ft^,  this  volume  of  air,  which  is  to  be  collected 
during  inflation,  is  calculated  from  the  steady-state  inflated 
shape  geometry  of  the  particular  parachute  type.  The  nomenclature 
is  described  in  Figure  24,p.A-14.  When  DQ  or  Dp  is  known,  a  can  be 
calculated  from  data  in  Table  I  and  Table  II,  p.  A-15,  for  various 
parachute  types  and  number  of  gores.  Then  the  geometric  volume  V0 
can  be  calculated  by  Equation  (31),  p.  A-14,  with  appropriate 
values  of  b/a  and  b' /a  from  the  tables. 

6.  AMq,  ft^,.  steady-state  canopy  mouth  area 


(B-1 ) 
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where  N/a,  b/a,  and  b'/a  are  available  from  Tables  I  and  II  for  the 
particular  type  of  parachute  and  number  of  gores. 

9  TT  2 

7.  Aso,  ft^,  canopy  surface  area  =  -^-Dq 

8.  C  ,  pressure  coefficient,  see  Figure  18,  p.  A-12.  A 
constant  Cp  =  1.7  for  all  altitudes  seems  to  yield  acceptable 
results . 


9.  Constants  k  and  n  are  derived  from  measurements  of  the  air 
flow  through  the  cloth.  Only  k  is  needed  for  Equation  (14),  but  n 
is  also  required  for  Equation  (13) .  These  parameters  can  be  deter¬ 
mined  for  any  cloth  using  the  technique  described  beginning  on 
p.  A-12.  The  two-point  method  is  adequate  if  the  AP  across  the 
cloth  is  in  the  range  of  AP  for  actual  operation.  Check-points  of 
cloth  permeability  can  be  measured  and  compared  to  calculated 
values  to  verify  agreement.  If  the  data  are  to  be  extrapolated  to 
operational  AP ' s  greater  than  measured,  a  better  method  of  deter¬ 
mining  k  and  n  from  the  test  data  would  be  a  least  squares  fit 
through  many  data  points.  This  way  errors  due  to  reading  either 
of  the  two  points  are  minimized. 

11.  Step  1 

Calculate  the  reference  time  tQ  by  use  of  Equations  (13)  or 
(14),  p.  A-7.  If  the  deployment  altitude  is  50,000  feet  or  higher. 
Equation  (14)  is  preferred  due  to  its  simplicity.  For  altitudes 
from  sea  level  to  50,000  feet.  Equation  (13)  is  preferred.  Figure 

12,  p.  A-8,  shows  the  effect  of  altitude  on  tQ  and  can  be  taken  as 

a  guide  for  the  user  to  decide  whether  to  use  Equation  (13)  or  (14). 
One  should  keep  in  mind  that  the  opening  shock  force  can  be  a 
strong  function  of  inflation  time,  so  be  as  realistic  as  possible. 

If  Equation  (13)  is  elected,  the  method  in  use  at  the  NSWC/WO  is 
to  program  Equation  (13)  to  compute  the  parachute  volume,  Vq,  for 
an  assumed  value  of  tQ.  Equation  (14)  ,  because  of  its  simplicity, 
can  be  used  for  a  first  estimate  of  tQ  at  all  altitudes.  The  com¬ 
puted  canopy  volume  is  then  compared  to  the  canopy  volume  calculated 
from  the  geometry  of  the  parachute  as  per  Equation  (31),  p.  A-14. 

If  the  volume  computed  from  the  mass  flow  is  within  the  volume 
computed  from  the  geometry  within  plus  or  minus  a  specified  delta 
volume,  the  time  tQ  is  printed  out.  If  not  within  the  specified 
limits,  tQ  is  adjusted,  and  a  new  volume  calculated.  For  a  35-foot 
D0 ,  T-10  type  canopy,  I  use  plus  or  minus  10  cubic  feet  in  the 
volume  comparison.  The  limit  would  be  reduced  for  a  parachute  of 
smaller  DQ. 


If 

VG 

calculated 

=  v0 

geometry 

± 

10, 

then 

print  answer. 

If 

v0 

calculated 

*  vo 

geometry 

+ 

10, 

then 

correct  tQ  as  follows 

■Vo  geometry 

Xo  “  *o  Y0  calculated  (B-2) 
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The  new  value  of  tD  is  substituted  in  the  'do  loop"  and  the 
volume  recomputed.  This  calculation  continues  until  the  required 
volume  is  within  the  specified  limits. 

III.  Calculate  tG  corrected  for  initial  area.  The  tQ  of  Section 
II  assumes  that  the  parachute  inflated  from  a  zero  initial  area. 

If  this  is  a  reasonable  assumption  for  the  particular  system  under 
study,  then  the  mass  ratio  can  be  determined  from  Equation  (6), 

p.  A-4.  For  ri  =  0  if  the  value  of  M  <  0.19,  then  a  finite  state 
of  deployment  exists,  and  the  time  ratio  of  occurrence  and  the 
maximum  shock  factor  can  be  determined  from  Equations  (9)  and  (10), 
respectively,  on  p.  A-5.  If  n  ^  0,  then  t-he  limiting  mass  ratio 
for  finite  operations  will  rise  slightly  as  described  in  Appendix 
C.  Figures  C-l  and  C-2  illustrate  the  effects  of  initial  area  on 
limiting  mass  ratios  and  shock  factors  respectively.  If  the  mass 
ratio  is  greater  th  a  n  the  limiting  mass  ratio  (Ml)  ,  then  the 
maximum  shock  force  occurs  at  a  time  greater  than  tQ  and  the 
elasticity  of  the  materials  must  be  considered  (see  Section  VI). 

If  n  /  0 ,  then  the  reference  time,  tQ ,  will  be  reduced,  and 
the  mass  ratio  will  rise  due  to  partial  inflation  at  the  line 
stretch.  Figures  9  and  10,  p.  A-6,  illustrate  the  effects  of 
initial  area  on  the  velocities  and  shock  factor  during  the 
"unfolding"  inflation.  Equation  (15),  p.  A-9,  can  be  used  to 
correct  tQ  calculated  for  the  cases  where  n  =  Ai/Ac.  If  the  initial 
value  of  drag  area  is  known,  Equation  (16),  p.  A-9,  can  be  used  to 
correct  tD  and  rechecked  for  limiting  mass  ratios  versus  n  in 
Appendix  C. 

IV.  Opening  shock  calculations  in  the  elastic  phase  of  inflation. 

It  has  been  considered  that  from  time  t  =  0  to  t  =  tQ  the  para¬ 
chute  has  been  inelastic.  At  the  time  t  =  tQ  the  applied  aero¬ 
dynamic  load  causes  the  materials  to  stretch  and  the  parachute 
canopy  increases  in  size.  The  increased  size  results  in  an 
increase  in  load,  which  causes  further  growth,  etc.  This  sequence 
of  events  continues  until  the  applied  forces  have  been  balanced  by 
the  strength  of  materials.  The  designer  must  insure  that  the 
constructed  strength  of  the  materials  is  sufficient  to  resist  the 
applied  loads  for  the  material  elongation  expected.  Use  of 
materials  of  low  elongation  should  result  in  lower  opening  shock 
forces  as  CDSmax  is  reduced. 

When  the  mass  ratio  of  the  system  is  greater  than  the  limiting 
mass  ratio,  the  elasticity  of  the  materials  and  material  strength 
determine  the  maximum  opening  shock  force.  The  maximum  elongation 
emax  an<^  the  ultimate  strength  of  the  materials  are  known  from 
tests  or  specifications.  The  technique  begins  on  p.  A-9. 

At  the  time  t  =  tD,  calculate  the  following  quantities  for  the 
particular  values  of  M  and  n . 

a.  vq/vs  from  Equation  (18),  p.  A-9. 
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b.  XQ  from  Equation  (22),  p.  A-10. 

c.  cQ  from  Equation  (23),  p.  A-10. 

d.  Determine  CDSmax/CpSQ  from  Figure  15,  p.  A-10. 

e.  Calculate  the  inflation  time  ratio  tf/tQ  from  Equation 
(24) ,  p.  A-10. 

f.  Calculate  the  maximum  shock  factor  from  Equation  (19), 

p .  A- 9 . 

g.  Calculate  the  opening  shock  force  Fm«x.  =  XiFs  where 


h.  Calculate  filling  time,  tf(sec) 


V.  In  order  to  simplify  the  required  effort,  the  work  sheets  of  Table 
B-l  are  included  on  pages  B-5  through  B-9  to  aid  the  engineer  in 
systematizing  the  analysis.  The  work  sheets  should  be  reproduced  to 
provide  additional  copies. 
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Table  B-l .  Opening  Shock  Force 


CALCULATION  WORK  SHEETS 


1.  Parachute  type  - 

2.  System  parameters 

a.  System  weight,  W  (lb) 

2 

b.  Gravity,  g  (ft/sec  ) 

c.  Deployment  altitude  (ft) 

d.  Deployment  air  density,  p  (slugs/ft^) 

e.  Velocity  at  line  stretch,  Vg  (fps) 

f.  Steady  state  canopy  data 


2a 


(1) 

Diameter,  DQ  (ft) 

(2) 

Inflated  diameter,  2a 

(ft); 

(3) 

Surface  area,  Sq  (ft^) 

’  4 

(4) 

Drag  area,  CDSQ  (ft2) ; 

CD 

(5) 

*  2 

Mouth  area,  A^q  (ft  ) 

MO 


.-2  [  /N/"a  -  b/aV 

*■  rl~) 


(6)  Volume,  V0”  (ft^) 

[*♦*] 


g.  Cloth  data 

(1)  k  )  Calculate  using  technique  beginning  on 


Note:  Permeability  is  usually 

measured  as  ft^/ft^/min. 

For  these  calculations 
permeability  must  be 
expressed  as  ft^/ft^/sec 


(2)  n  )  p.  A-12. 


k 

n 


*  Data  for  these  calculations  are  listed  in  Tables  1 
and  2,  p.  A-15. 
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Table  B-l  .  Opening  Shock  Force  N'J'k  N& 

(cont'd)  'A 

(3)  emax»  determine  maximum  elongations  from  ' 

pull  test  data  of  joints,  seams,  lines,  etc.  Use 
minimum  Smax  determined  from  tests. 


(4)  Cp;  pressure  coefficient 
h.  Steady  state  drag,  Fs  (lb),  Fs  -  CDSo 


F. 


i.  Parachute  constructed  strength,  Fc  (lb);  deter¬ 
mined  from  data  on  efficiency  of  seams,  joints,  lines. 
Constructed  strength  is  the  minimum  load  required  to 
fail  a  member  times  the  number  of  members. 


3.  Force  calculations 

a.  Calculate  tQ  for  n  =  0;  eq.  14,  p.  A-7 « 


paYo 


14W 


P9VSCDSo 


CDSo 


2W  L  AMO  ”  Aso  k  /CpP 


Check  Figure  13,  p.  A-8,  for  advisability  of 
using  eq.  13,  p.  A-7. 

b.  If  n  =  0,  proceed  with  steps  c  through  e. 

If  ri  ^  0,  go  to  step  f. 

c.  Mass  ratio,  M;  eq.  6,  p.  A-4 

2W 

M  =  - 

P9Vst0CDSo 


M 


d.  If  M  <  4/21  for  n  =  0,  then  finite  mass 
deployment  is  indicated. 

(1)  Time  ratio  at  x^  max;  eq.  9,  p.  A-5 


1 


(2)  Max  shock  factor,  x^;  eq. ,  10,  p.  A-5 
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Table  B-l .  Opening  Shock  Force 

(Cont'd) 

(3)  Max  shock  force,  F  (lb) 

uld  A 


F  =  X. 

max.  i  , 


e.  If  M  >  4/21;  then  intermediate  mass  or  infinite 
mass  deployment  is  indicated  and  the  elasticity  of 
materials  is  involved.  Calculate  the  trajectory  con¬ 
ditions  at  time  t  =  t^. 


(1)  Velocity  ratio  @  t 


p.  A- 10 


19,  p.  A- 9 


tQ  for  n  =  0 


V 


1  + 


7M 


(2)  Shock  factor  XQ  @  t  =  tQ  for  r] 


V 


=  0 


['**] 


2  l  V0 


(3)  Initial  elongation,  eQ;  eq .  23,  p.  A-10 


X  F- 

e.  =  — 2-S-  e„ 


cns 

(4)  Determine  - from  Figure  15,  p.  A-10 

CDSo 

(5)  Calculate  inflation  time  ratio,  tf  ;  eq.  24, 


**DSmax 


(6)  Calculate  maximum  shock  factor,  max;  eq 


O' 


X: 


_vo 


7  — ,  — 1 2 

— 1 


(7)  Calculate  maximum  shock  force,  Fmax  (lb). 


Fmax=  Xj  FS 
max 


V 


cDSmax 

cDso 


tf 

to 


X  : 
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Table  B-l .  Opening  Shock  Force 

(cont ' d) 


(8)  Inflation  time,  sec  =  tf  =  t. 


■d) 


f.  If  r|  ^  0,  correct  tQ  for  initial  area  effects; 

eq.  16,  p.  A- 9  ("  /C^S.V76] 

to  *  r*  FT  )  to  ca,cu,ated 

g.  Mass  Ratio,  M,  eq.  6,  p.  A-4  \D  o/  J 


M 


2W 


h.  Calculate  limiting  mass  ratio,  M^- 

i  r  9  2  3  1 1 

ml  =  3iT^r[uT?  *ui+7j 


If  M  <  Mt  ,  finite  mass  deployment  is  indicated  and  x, 

u  i  max 

can  be  determined  by  eq .  8,  p.  A-5  by  assuming  values 

of  t/tQ  and  plotting  the  data  using  the  methods  of 

Appendix  C. 

i.  If  M  >  M^,  then  intermediate  mass  or 
infinite  mass  deployment  is  indicated  and  the 
elasticity  of  materials  is  involved.  Calculate 
the  trajectory  conditions  at  time  t  =  tG. 


p.  A- 9 


(1)  Velocity  ratio  @  t  =  tQ  for  n  t  0;  eq.  18, 


V. 


22,  p.  A- 10 


VS 


(2)  Shock  factor  X  @  t  =  t  for  n  4  0;  eq. 


rV 


—  VVs 


(3)  Initial  elongation,  eo;  eq.  23,  p.  A-10 


XoFS 


C„S 


(4)  Determine  ..D  max  fr0m  Figure  15,  p.  A-10 


CDSo 


M 


M, 


V„ 


CDsi 


max 


CDS0 


Sac. 


Sec. 
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Table  B-l.  Opening  Shock  Force  {Con 


(5)  Calculate  inflation  time  ratio,  —  ;  eq.  24, 
p.  A- 10  -  to 


CDSm.x 


(6)  Calculate  maximum  shock  factor,  max; 
eq.  19,  p.  A- 9  /t.\6 


nax  r  V, 


R-MO-HI 


(7)  Calculate  maximum  shock  force,  Fmax(lb) 


Fm.x  =  Xi  FS 
max 


(8)  Calculate  inflation  time,  tf(sec) 

'--m 
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